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I review recent developments in lattice QCD. I first give an overview of its formalism, 
and then discuss lattice discretizations of fermions. We then turn to a description of the 
quenched approximation and why it is disappearing as a vehicle for QCD phenomenology. 
I describe recent claims for progress in simulations which include dynamical fermions 
and the interesting theoretical problems they raise. I conclude with brief descriptions of 
the calculations of matrix elements in heavy flavor systems and for kaons. 
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1. Introduction 

Lattice regulated quantum field theory supplemented by numerical simulations is a 
major source of information about Quantum Chromodynamics (QCD), the theory 
of the strong interactions. This article is an overview of lattice QCD, with a target 
audience of particle physicists who do not do lattice calculations themselves, but 
might be "users" of them. Since lattice QCD is a mature field El most work in lat- 
tice QCD represents evolutionary progress, not revolutionary breakthroughs. But 
lattice QCD is presently going through a transition phase, as a long-used approx- 
imation to QCD (the quenched approximation) becomes increasingly problematic, 
and numerical studies of QCD which include the effects of dynamical quarks are 
beginning to produce interesting numbers - and their own interesting questions of 
principle. 

Most lattice reviews -like most lattice calculations - are fixated on numbers, 
the values of reduced matrix elements computed at physical values of quark masses. 
Numbers are important, and the lattice can be the best source of some specific 
(model-independent) predictions of QCD. However, numbers evolve and become 
obsolete. I think it will be more interesting to take a more impressionistic view 
of lattice QCD, and concentrate on how the calculations are done. There are a 
lot of ways that they can fail, and phcnomcnologists who are going to use lattice 
determinations of reduced matrix elements in their own work ought to be aware of 
the potential weaknesses in them. 

I will begin with a quick overview of lattice methodology: how to go from a con- 
tinuum action to a lattice calculation, and then get back to the continuum. Next, I 
will describe methods for putting fermions on the lattice. A lot has happened here 
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since the last textbooks were written. I will then describe the quenched approxi- 
mation (which is about to disappear as a source of quantitative lattice predictions) 
and tell you about recent work with dynamical fermions, which is certainly in a 
very peculiar state. I will then give an overview of two kinds of standard model 
tests, first heavy flavor physics, and then weak matrix elements of kaons. I will not 
discuss QCD thermodynamics, lattice-motivated scenarios of confinement or chiral 
symmetry breaking, or iV c ^ 3. 

2. A superficial overview of lattice methodology 
2.1. Lattice Variables and Actions 

In order to perform calculations in quantum field theory it is necessary to control 
their ultraviolet divergences. The lattice is a space-time cutoff which eliminates 
all degrees of freedom from distances shorter than the lattice spacing a. As with 
any regulator, it must be removed after renormalization. Contact with experiment 
only exists in the continuum limit, when the lattice spacing is taken to zero. Other 
regularization schemes are tied closely to perturbative expansions: one calculates 
a process to some order in a coupling constant; divergences are removed order by 
order in perturbation theory In contrast, the lattice is a non-perturbative cutoff. 
Before a calculation begins, all wavelengths less than a lattice spacing are removed. 

All regulators have a price. On the lattice we sacrifice all continuous space-time 
symmetries but preserve all internal symmetries, including local gauge invariance. 
This preservation is important for non-perturbative physics. For example, gauge 
invariance is a property of the continuum theory which is non-perturbative, so 
maintaining it as we pass to the lattice means that all of its consequences (including 
current conservation and renormalizability) will be preserved. Of course, we want 
to make predictions for the real world, when all symmetries are present. It can be 
very expensive to recover the symmetries we lost by going on the lattice. 

It is straightforward to construct the lattice version of a field theory. One just 
replaces the space-time coordinate x^ by a set of integers (x M = an M , where 
a is the lattice spacing). Field variables (f>(x) are defined on sites <p(x n ) = (j> n , 
The action, an integral over the Lagrangian, is replaced by a sum over sites, and 
derivatives in the Lagrange density are replaced by finite differences. 

S = f d 4 x£(t(x)) - a 4 ]T £(M (1) 

and the generating functional for Euclidean Green's functions is replaced by an 
ordinary integral over the lattice fields 

Z= /(Q^™)e s . (2) 

J n 

Gauge fields carry a space-time index /i in addition to an internal symmetry index 
a (A^(x)) and are associated with a path in space x fl (s): a particle traversing a 
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contour s in space picks up a phase factor 

i\> -> P(exp ig J dx^A^tp = U(s)tp(x). (3) 

P is a path-ordering factor analogous to the time-ordering operator in ordinary 
quantum mechanics. Under a gauge transformation g, U(s) is rotated at each end: 

U(s) -> n^MsMMO)- 1 ). (4) 

These considerations led Wilson to formulate gauge fields on a space-time lattice, 
in terms of a set of fundamental variables which are elements of the gauge group 
G (I'll specialize to SU (N) for TV colors) living on the links of a four-dimensional 
lattice, connecting neighboring sites x and x + a/i: U^{x), 

U^x) = exp^aTM^x)) (5) 

(g is the coupling, the vector potential, and T a is a group generator). 
Under a gauge transformation link variables transform as 

U^x) -> V(x)U„(x)V(x + A) f (6) 

and site variables transform as 

ip(x) -> V{x)ip{x). (7) 

We are usually only interested in gauge invariant observables. These will be either 
matter fields connected by oriented "strings" of U's 

ip(x 1 )U fl (xi)U^(x 1 +p,)...ip(x 2 ) (8) 

or closed oriented loops of U's 

Tr . . . U^x)U^{x + A) ■ . . Tr . . . U„(x)V^{x + ft)V(x + ftU^x + A) ■ • ■ ■ (9) 

An action is specified by recalling that the classical Yang-Mills action involves 
the curl of A^, F^ v . Thus a lattice action ought to involve a product of U^s around 
some closed contour. Gauge invariance will automatically be satisfied for actions 
built of powers of traces of U's around arbitrary closed loops, with arbitrary coupling 
constants. If we assume that the gauge fields are smooth, we can expand the link 
variables in a power series in gaA'^s. For almost any closed loop, the leading term in 
the expansion will be proportional to F^ v . All lattice actions are just bare actions 
characterized by many bare parameters (coefficients of loops). In the continuum 
(scaling) limit all these actions are presumed to lie the same universality class, 
which is (presumably) the same universality class as QCD with any regularization 
scheme, and there will be cutoff-independent predictions from any lattice action 
which are simply predictions of QCD. 

The simplest contour has a perimeter of four links. The "plaquette action" or 
"Wilson action" (after its inventor) is defined as 

S = \ E E Re Tr i 1 - U MMn + + 0)Ul{n)). (10) 
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The bare lattice coupling, whose associated cutoff is a, is g 2 . The lattice parameter 
(3 — 2N/g 2 is often written instead of g 2 = 47ra s . 

In the strong coupling limit, the lattice regularized version of a gauge theory 
with any internal symmetry group automatically confined^. The interesting question 
is to understand what happens as the bare coupling became weaker at fixed lattice 
spacing (or for asymptotically free theories, what happens as the lattice spacing is 
taken to zero). In the strong coupling limit, chiral symmetry is also spontaneously 
brokerM While it is still not proved, all the evidence we have is that if the internal 
symmetry group is SU(N), confinement and chiral symmetry persist as the lattice 
spacing is taken away. 

2.2. Numerical Simulations 

Lattice QCD has survived because it is a framework for doing "exact" (direct from 
a cutoff Lagrangian) calculations of QCD, which is simple enough in principle that 
it can be taught to a computer. In a lattice calculation, like any other calculation 
in quantum field theory, we compute an expectation value of any observable T as 
an average over a ensemble of field configurations: 



The average is done by numerical simulation: we construct an ensemble of states 
(collection of field variables), where the probability of finding a particular con- 
figuration in the ensemble is given by Boltzmann weighting (i. e. proportional to 
exp(— S). Then the expectation value of any observable T is given simply by an 
average over the ensemble: 



As the number of measurements N becomes large the quantity T will become a 
Gaussian distribution about a mean value, our desired expectation value. 

The idea of essentially all simulation algorithm^ to construct a new config- 
uration of field variables from an old one. One begins with some initial field con- 
figuration and monitors observables while the algorithm steps along. After some 
number of steps, the value of observables will appear to become independent of the 
starting configuration. At that point the system is said to be "in equilibrium" and 
Ea. 1121 can be used to make measurements. 

Dynamical fermions are a complication for QCD. The fermion path integral is 
not a number and a computer can't simulate fermions directly. For n/ degenerate 
fermion flavors the generating functional for Green's functions is 




(11) 




(12) 



i=l 



Z= [dU}[di/j][dtl;] exp{~(3S G (U) -^2ipM(U)ip) 



(13) 
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One formally integrates out the fermions to give a pure gauge action where the 
probability measure includes a nonlocal interaction among the C/'s: 

Z= f[dU] (detAf (U)) n f exp(-pS(U)) (14) 



[dU] enp(-PS(U) + n/Trm(M(£/))). (15) 

Generating configurations of the U's involves computing how the action changes 
when the set of C/'s are varied. The presence of the determinant makes this problem 
very difficult. For a pure gauge theory, changing a variable at one location only 
affects the action at sites "near" the variable, so the attempt to update one link 
variable on the lattice involves a computational effort independent of the lattice 
volume (said differently, the cost of generating a new configuration scales with the 
lattice volume). However, the determinant is nonlocal, and so in principle updating 
one gauge variable would involve an amount of work proportional to the lattice 
volume. 

Typically, this involves inverting the fermion matrix M(U) (because the change 
in logM is d log M/dU = M~ 1 dM/dU). This is the major computational problem 
dynamical fermion simulations face. M has eigenvalues with a very large range- 
from 2ir down to m q a. The "conditioning number" - the ratio of largest to smallest 
eigenvalue of the matrix - determines the convergence rate of iterative methods 
which invert it. In the physically interesting limit of small m q the conditioning 
number diverges and the matrix is said to become "ill-conditioned." The matrix 
becomes difficult (impossible) to invert. At present it is necessary to compute at 
unphysically heavy values of the quark mass and to extrapolate to m q — 0. (The 
standard inversion technique today is one of the variants of the conjugate gradient 
algorithm. ) 

The tremendous expense of the determinant is responsible for one of the stan- 
dard lattice approximations, the "quenched" approximation. In this approximation 
the back-reaction of the fermions on the gauge fields is neglected, by setting nj — 
in Eq. 1131 Valence quarks, or quarks which appear in observables, are kept, but no 
sea quarks. 

All large scale dynamical fermion simulations today generate configurations us- 
ing some variation of the microcanonical ensemble. That is, they introduce momen- 
tum variables P conjugate to the C/'s and integrate Hamilton's equations through 
a simulation time t 

U = iPU (16) 

The integration is done numerically by introducing a time step At. The momenta 
are repeatedly refreshed by bringing them in contact with a heat bath and the 
method is thus called Refreshed or Hybrid Molecular Dynamics^. The reason for 
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the use of these small time step algorithms is that for any change in any of the U's, 
M^ 1 must be recomputed. When Eon. 1171 is integrated all of the C/'s in the lattice 
are updated simultaneously, and only one matrix inversion is needed per change of 
all the bosonic variables. 

For special values of n/ the equations of motion can be derived from a local 
Hamiltonian and in that case At systematics in the integration can be removed by 
an extra Metropolis accept/reject step. This method is called Hybrid Monte Carlo™ 
There are also many variations on these schemed. 

Recently, there has been a lot of exploration of alternative, large-step algorithms 
for dynamical fermionJ^. While these look promising, none has been used as much 
as the molecular dynamics methods. 



2.3. Spectroscopy Calculations 

Masses of hadrons are computed in lattice simulations from the asymptotic behav- 
ior of Euclidean-time correlation functions. A typical (diagonal) correlator can be 
written as 

K{x) = (0|O(x)O(0)|0). (18) 
This correlator has its Euclidean-space analog 

U(q) = J d 4 xcxp{iqx)K(x) (19) 

and if we assume that n(q) is dominated by a sum of resonances, we expect 

' q z + mz 

where r n = |(0|O|n)| 2 . We can invert the Fourier transform to obtain a prediction 
for K{x) (or some re- weighted analog). For example, lattice people often sum K{x) 
over a three-dimensional volume of points (x = {t, r}, sum r) to produce 

C{t)=Y J K{x) 



dqo 

e-™" 4 (21) 



\ " r " r ~rn n t 



^ 2m„ 

n 

At large separation the correlation function is approximately 

C{t) = J^e— * (22) 

where mi is the mass of the lightest state which the operator O can create from the 
vacuum. Fig^shows an example of a lattice correlator. Obviously, the lightest state 
in a channel is the easiest one to study. Excited states are hard to work with (their 
signal goes under the ground state's.) If the operator has vacuum quantum numbers 
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Fig. 1. A "better than typical" lattice correlator and its fit. Periodic boundary conditions convert 
the exponential decay into a hyperbolic cosine. 

(an operator coupling to the scalar glucball is an example), the contribution for the 
"real" lightest state will disappear exponentially into the constant background. In 
that case it is also hard to measure a mass. 

Most of the observables we are interested in will involve valence fermions. To 
compute the mass of a meson we might take 

C{t)=^2{J{x,t)J{0,0)) (23) 

where 

J(x, t) = tp(x, t)Tip(x, t) (24) 

and r is a Dirac matrix. The intermediate states \n) which saturate C(x, t) are 
the hadrons which the current J can create from the vacuum: the pion, for a pseu- 
doscalar current, the rho, for a vector current, and so on. Now wc write out the 
correlator in terms of fermion fields 

c(t) = 5>|$« (M)r«^(M)$(o,o)r fcI ^f (o,o)|o> (25) 

X 

with a Roman index for spin and a Greek index for color. We contract creation and 
annihilation operators into quark propagators 

(0|T(V>, t)$(0, 0))|0> = Gf{x, t; 0, 0) (26) 

so 

C(t) = TrG(x, t; 0, 0)rG(0, 0; x, t)T (27) 

X 
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where the trace runs over spin and color indices. Baryons are constructed similarly. 

In practice, lattice people do not use such simple operators for interpolating 
fields. Instead, they take more complicated operators which model the wave function 
of the meson. One way to do this would be to gauge fix the lattice to some smooth 
gauge, like Coulomb gauge, and take O(xo,t) = J2 X y ip(x,t)r?p(y,t)§(x — x ,y — 
xq), where $ is our guess for the wave function (I like Gaussians). These guesses 
often have parameters which can be "tuned" to enhance the signal. It's important 
to note that the choice of trial function should not affect the measured mass, it just 
alters how quickly the leading term emerges from the sum. 

2.4. Getting rid of the lattice 

The lattice spacing a is unphysical: it was an UV cutoff we introduced by hand 
to regularize the theory. The couplings {g}we typed into the computer were bare 
couplings which define the cutoff theory at cutoff scale a. When we take a to zero 
we must also specify how the couplings behave. The proper continuum limit comes 
when we take a to zero holding physical quantities fixed, not when we take a to 
zero at fixed lattice couplings {<?}. (It's conventional to think of a as a function of 
the bare couplings {<?}, so people actually think about tuning the {g}'s to push a to 
zero. A lattice QCD action will have one marginally relevant coupling-the lattice 
analog of the gauge coupling - plus a set of coefficients of irrelevant operators. 
When they are present, quark masses are relevant operators, and they must also be 
tuned as the lattice spacing is varied. These quantities characterize the particular 
lattice discretization.) 

On the lattice, if all quark masses are set to zero, the only dimensionful pa- 
rameter is the lattice spacing, so all masses scale like I/a. Lattice Monte Carlo 
predictions are of dimcnsionlcss ratios of dimensionful quantities. One can deter- 
mine the lattice spacing by fixing one mass from experiment, and then one can go 
on to predict any other dimensionful quantity. However, at nonzero lattice spacing, 
all predictions will depend on the lattice spacing. Imagine computing some masses 
at several values of the lattice spacing. (Pick several values of the bare parame- 
ters and calculate masses for each set of couplings.) If the lattice spacing is small 
enough, the typical behavior of a ratio will look like 

(omi(a))/(am2(a)) = mi(0)/m2(0) + 0{m\a ) + 0(( mi a) 2 ) + . . . (28) 

(modulo powers of log (mi a)). The leading term does not depend on the value of the 
UV cutoff. That is our cutoff-independent prediction. Everything else is an artifact 
of the calculation. We say that a calculation "scales" if the a— dependent terms in 
Eq. [28] are zero or small enough that one can extrapolate to a = 0, and generically 
refer to all the a— dependent terms as "scale violations." Clearly our engineering 
goal is to design our lattice simulation to minimize scale violations. 

We can imagine expressing each dimensionless combination am(a) as some func- 
tion of the bare coupling(s) {g(a)}, am — f({g(a)}). As a — + we must tune the 
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set of couplings {g(a)} so 

lim -f({g(a)}) — * constant. (29) 
a^o a 

From the point of view of the lattice theory, we must tune {g} so that correlation 
lengths 1/ma diverge. This will occur only at the locations of second (or higher) 
order phase transitions. In QCD the fixed point is g c = so we must tune the 
coupling to vanish as a goes to zero. 

One needs to set the scale by taking one experimental number as input. A 
complication that you may not have thought of is that the theory we simulate on 
the computer might be different from the real world, so its spectrum would be 
different. For example, the quenched approximation, or for that matter QCD with 
two flavors of degenerate quarks, almost certainly does not have the same spectrum 
as QCD with six flavors of dynamical quarks with their appropriate masses. Using 
one mass to set the scale from one of these approximations to the real world might 
not give a prediction for another quantity which agrees with experiment. 

2.5. Lattice Error Bars 

Numbers presented from lattice simulations come with uncertainties. Phenome- 
nologists ought to read carefully the parts of the papers which describe the error 
analysis, because there are many parts to a lattice number's uncertainty, all differ- 
ent. Some of the uncertainty is statistical: The sample of lattices is finite. Typically, 
the quoted statistical uncertainty includes uncertainty from a fit: it is rare that a 
simulation measures one global quantity which is the desired observable. Usually 
one has to take a lattice correlator, assume it has some functional form, and de- 
termine the parameters which characterize the shape of the curve. The fit function 
fea. !22l is an example) may be only asymptotic, and one has to figure out what part 
of the data is described by one's function. A complication which enters in at this 
point is that different quantities measured on the same set of lattices are typically 
highly correlated. These correlations have to be taken into account in the fit. 

Dimensionful quantities require (at least) two measurements: an additional 
quantity is needed to set the scale. Again, different quantities may be easier or 
harder to measure, so their lattice errors will vary. Clearly, people like to choose 
quantities which are as insensitive as possible to lattice artifacts, interpolation in 
mass, or other analysis issues. For example, one might want to use a hadron mass, 
or a decay constant, or a parameter associated with the heavy quark potential, to 
fix a. Many people use the rho mass, but perhaps this is not a good quantity: in 
the real world, the rho is unstable (and broad), so what will happen in a simu- 
lations? The S-P mass splitting in mesons is known to be relatively independent 
of the quark mass, but for light hadrons, the P-wave signal is noisy. Using heavy 
quark properties to fix the lattice spacing usually requires a different heavy quark 
action, with a separate simulation. There are endless arguments in the literature... 
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This is not the end of the analysis. It may be necessary to extrapolate or in- 
terpolate one's results to a particular value of the quark mass. One again needs 
a functional form. The fit or extrapolation has a statistical uncertainty, but now 
systematics begin to creep in: The functional dependence comes from some theory. 
Is the theory correct? Can it be applied to the whole range of quark masses where 
the lattice data exists? 

Additional systematics arise in quenched approximation. Quenched QCD is 
not real- world QCD, and their spectra are presumably different. People publish- 
ing quenched numbers sometimes attempt to quote a systematic uncertainty from 
the quenched approximation. They might try to do this by using several different 
observables to set the lattice spacing, or to fix the quark mass, and seeing how 
the final result changes for different observables. The problem with this analysis 
is that it assumes that if the parameter values were chosen in the same way in 
quenched and full QCD (take the strange quark mass from the K* /<f> mass ratio, 
for example), the desired matrix element would also be the same in quenched and 
full QCD. There is no reason for that to be true. In full QCD there is supposed to 
be one unique value of the lattice spacing, so differences from different observables 
reflect on the quality of the simulation (assuming that QCD is in fact the correct 
description of Nature) . 

Finally, lattice quantities which are not spectral have scheme dependence. It is 
necessary to convert the lattice number to a number in some continuum regular- 
ization scheme. If this scheme matching is done nonperturbatively, the uncertainty 
in doing this is likely to be mostly statistical: the conversion factor comes from 
its own simulation. When the matching factor is computed in perturbation theory, 
most people will make "reasonable" choices for the way the calculation is organized 
(picking the scale at which the strong interaction coupling constant is determined, 
for example), vary their choices in some "reasonable" way, and attempt to assign 
an uncertainty based on the variation they see. 

It is relatively easy to do simulations of gauge theories for any compact internal 
symmetry group, or any other bosonic system. Fermionic systems interacting with 
gauge or matter fields are feasible without enormous resources for fermion masses 
down to the strange quark mass, but become quite expensive below that value. 
Simulations with massless or nearly massless fermions remain at the costly frontier. 
Of course, large N c simulations scale roughly like iV 3 from the multiplication of the 
link matrices. In spectroscopy, properties of flavor non-singlet hadrons are easier to 
compute than those of flavor singlet ones (disconnected diagrams are noisy), and 
almost any calculation can be designed to scale linearly with the number of quark 
propagators which need to be strung together. 

2.6. Improvement: Why and How 

Today's QCD simulations range from 16 3 x 32 to 32 3 x 100 points and run from 
hundreds (quenched) to thousands (full QCD) of hours on the fastest supercom- 
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puters in the world. The cost of a Monte Carlo simulation in a box of physical size 
L with lattice spacing a and quark mass m q scales roughly as 

(-) 4 (-r 2 (-r 3 (30) 

a' a' ra q ' 

where the first term in parentheses just counts the number of sites, the second term 
gives the cost of "critical slowing down" -the extent to which successive configu- 
rations are correlated, and the third term gives the cost of inverting the fermion 
propagator, plus critical slowing down from the nearly massless pions. Thus it is 
worthwhile to think about how to do the discretization, to maximize the value of 
the lattice spacing at which useful (scaling) simulations can be done. 

Remember, the lattice action is just a bare action defined with a cutoff. No 
lattice discretization is any better or worse (in principle) than any other. Any 
bare action which is in the same universality class as QCD will produce universal 
numbers in the scaling limit. However, by clever engineering, it might be possible 
to devise actions whose scaling behavior is better, and which can be used at bigger 
lattice spacing. 

To want to compute the value of some QCD observable via numerical simulation 
is to confront a daunting set of technical problems. The lattice volume must be 
large enough to contain the hadrons - and often, to allow them to be located far 
apart. High statistics are needed for reasonable accuracy. And the lattice spacing 
must be small enough to minimize discretization artifacts. These constraints push 
simulations onto large, fast, and often remote supercomputers. This introduces a 
new set of problems. These resources are expensive, so people form collaborations 
to share them. The fastest machines are the newest ones, which are often unstable, 
or hard to program, inefficient for all but restricted kinds of computations, or all 
of the above. Finally, because they are so big, lattice projects have a high profile. 
They cannot be allowed to fail. 

Thus, until the late '90's, most simulations used the simplest discretizations, 
or minimal modifications to them. The last few years have seen a slow change in 
this situation, basically driven by necessity: standard methods become increasingly 
dominated by artifacts as the light quark masses become smaller and smaller. Two 
examples illustrate this. 

First, groups doing dynamical fermion simulations with staggered fermions (see 
Sec. 3) want to do simulations with two nearly degenerate light quarks (the up and 
down quarks) and a heavier strange quark. With this kind of fermion, it is possible to 
do simulations with the strange quark mass at its physical value. However, staggered 
fermions suffer from flavor symmetry breaking (see Sec. 3.2), and at today's lattice 
spacings the standard formulation of staggered fermions produces a spectrum of 
mesons in which some of the non-strange pseudoscalars would be heavier than 
the strange ones. Then what is a pion, and what is a kaon? A more complicated 
discretization is needed to reduce the flavor symmetry breaking, so that the non- 
strange mesons are separated in mass from the strange ones. 
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A second example is the discovery of lattice discretizations with flavor symmetry 
and improved or exact chiral symmetry (domain wall and overlap fermions) . Their 
extra computational expense is repaid by their ability to cleanly address problems 
(mostly involving chiral symmetry) the standard methods cannot. 

The democratization of computing has helped. QCD is not yet possible on desk 
top machines, but a lot of interesting work can be done on clusters for a cost of a 
few tens of thousands of dollars. 

(Of course, most people will still minimize the amount of evolution of their codes 
away from the well-studied algorithms because the original requirements of large 
volume and high statistics have not gone away.) 

The simplest organizing principle for "improvement" is to use the canonical 
dimensionality of operators as a guide. Consider the gauge action as an example. 
If we perform a naive Taylor expansion of a lattice operator like the plaquette, we 
find that it can be written as 

1 - ^Rc TrU p i aq = roTrF^ + a 2 [n £ M „ Tr D^D^ + 

+0(a 4 ) (31) 

The expansion coefficients have a power series expansion in the bare coupling, 
rj = Aj+g 2 Bj + . . .. Other loops have a similar expansion, with different coefficients. 
The expectation value of any operator T computed using the plaquette action will 
have an expansion 

(T(a)> = (T(0)) + O(a) + 0(g 2 a) + ... (32) 

Now the idea is to take the lattice action to be a minimal subset of loops and 
systematically remove the a™ terms for physical observables order by order in n 
by taking the right linear combination of loops in the action, S = ^ • CjOj with 
Cj = c® + g 2 c^ + . . .. This method was developed by Symanzik and co-worker J^ | 10 | l 1 J 
in the mid-80's. 

Rcnormalization group ideas have also been used to motivate improvement 
programJi^l. 

If you think about it, all the operators which are being added to the minimal 
discretizations are irrelevant operators. So in principle, no improvement method 
is really "wrong," although it might happen that the cost of the improved action 
is greater t han the gain from simulation fidelity. The final arbiter is a simulation. 
An exampleE*! of a test of scale violations is shown in Fig. [21 The x axis is the 
lattice spacing, in units of a quantity r\, which is defined through the heavy quark 
potential: r\dV (r) / dr\ ri — 1.0, about 0.4 fm. The plotting symbols are for different 
kinds of discretizations. The flatter the curve, the smaller the scale violations. The 
data labeled by squares shows the best scaling. 
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Fig. 2. Lattice calculations of the (a) rho and (b) nucleon mass, interpolated to the point 
m^r\ = 0.778, as a function of lattice spacing, from Ref. 1131 



A recent innovation for improvement is the use of "fat links" as gauge connec- 
tions in fermion actions. A standard "thin link" action has its gauge connections 
built out of a single link variable, S ~ tp(x)U fJj (x)ip(x + ft) . Fat links replace th e sin - 
gle link by a average of paths. For example, a simple blocking ("APE blockin^Bl') 
is 

V fi (x)= (l-a)iyx) 

+«/6 E^(UAx)U^(x + v)U v (x + ^ 
+U v {x-v)^U^{x-i>)U v {x-v + (£)), (33) 

with possibly a further projection back to the gauge group. The parameter a can be 
tuned. There are of course many possibilities This year, the most popular choices 
for fattening include the Asqtad linlffSl (" a 2 TADpole improved") and the HYP 
(hypercubic) linlff^. 

For smooth fields the fat links have an expansion V^(ar) = 1 + iaB^(x) + . . . and 
the original thin links have an expansion U^(x) = 1 + iaA fi (x) + . . ., where 

B^x) = J2 h ^(v)Mx + y) (34) 

y," 

and the convolution function h^ v (y) obviously depends on the fattening. Smoothing 
is easiest to see in momentum space, w here the connection between the thin and 
fat link basically introduces a form factoi^^, B^(q) — J2„ h^ u (q)A v (q). 

A number of ideas motivate fat links. They all boil down to the simple idea 
that a hard lattice cutoff introduces more discretization artifacts than a smoothed 
cutoff. 

Clearly, the trade off is between smoothing the gauge field locally (at the cut- 
off scale) versus erasing physics at long distances. Fat links improve perturbation 
theory, reduce flavor symmetry breaking for staggered fermions, reduce chiral sym- 
metry breaking for non-chiral discretizations, and help in the implementation of 
overlap action*^. 
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3. Relativistic Fermions on the Lattice 

Finding a lattice discretization for light fermions involves yet another problem: 
doubling. Let's illustrate this with free field theory. The continuum free action is 

S = J d A x$(x) ltl d^(x) + m-4)(x)ip(x)]. (35) 

One obtains the so-called naive lattice formulation by replacing the derivatives by 
symmetric differences: we explicitly introduce the lattice spacing a in the denomi- 
nator and write 

gnaive = £ ^ A/1 ^ n + TO ^ (36) 

where the lattice derivative is 

The propagator is easy to construct: 

-G[p) = (17^ smp^a + ma) = — — — ■ — . (38) 

a sm P^ a ~i~ m a 

Now the lattice momentum p M ranges from — ir/a to n/a. A continuum fermion 
with its propagator (ijuPu + to) -1 has a large contribution at small p from four 
modes which are bundled together into a single Dirac spinor. The lattice propagator 
has these modes too, at p = (0, 0, 0, 0), but there are other degenerate ones, at p = 
(tt, 0, 0, 0), (0, 7r, 0, 0), ... (n, 7r, 7r, 7r). As a goes to zero, the propagator is dominated 
by the places where the denominator is small, and there are sixteen of these (64 
modes in all) in all the corners of the Brillouin zone. Thus our action is a model 
for sixteen light fermions, not one. This is the famous "doubling problem." 

The doubling problem is closely connected to the axial anomaly. Karsten and 
SmiPl showed by explicit calculations that the axial charges of the sixteen light 
fermions are paired and sum to zero. Many years ago Adler showed that it is 
not possible to find a continuum regulator which is gauge invariant for a theory 
with continuous chiral symmetry and a weak coupling limit with a perturbative 
expansiorP^, and the lat tice was not immune to this problem. 

Nielsen and Ninomiya^J codified this constraint in a famous "no-go" theorem. 
In detail, the theorem assumes 

• A quadratic fermion action ip(x)H(x — y)ip{y), where H is Hermitian, has 
a Fourier transform H{p) defined for all p in the Brillouin zone, and has 
a continuous first derivative everywhere in the Brillouin zone. H(p) should 
behave as for small p^ . 

• A local conserved charge Q defined as Q = J2 X 3o{ x )< where jo is a function 
of the field variables tp{y) where y is close to x. 

• Q is quantized. 
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The statement of the theorem is that, once these conditions hold, H(p) has an equal 
number of left handed and right handed fermions for each eigenvalue of Q. 

There is a folkloric version of the theorem which says that no lattice action 
can be undoubled, chiral, and have couplings which extend over a finite number of 
lattice spacings ("ultra-locality"). This is actually not what the theorem says. It is 
the quantization of the charge which governs whether the theorem is evaded or not. 

Ultra-locality is a historical engineering constraint on lattice action design. What 
is needed for a proper field theoretic description is locality, meaning that the range 
of the action is restricted to be on the order of the size of the spatial cutoff. It 
is believed that having lattice couplings which fall off exponentially with distance 
(measured in units of the lattice spacing), i.e. S — J2x r i>{x)C(r)i()(x + r) with 
C(r) ~ exp(— r/£), £ oc a, corresponds to a local action in the continuum limit, 
but that slower falloff (power law, for example) does not. Present day simulations 
have cutoffs of 0.2 to 0.05 fm, which is uncomfortably close to physical scales, and 
so people who work with non-ultra-local actions worry about the range of their 
actions, and try to tune them to maximize locality. 

At the end of their paper, Nielsen and Ninomiya discuss ways to evade the the- 
orem. They were not too optimistic about abandoning the quantization of charge, 
but that is how overlap and domain wall fermions achieve chirality. 

Having said that, it is useful to classify lattice actions into their folkloric cate- 
gories: ultra-local actions which are non-chiral and undoubled (Wilson fermions and 
their generalizations), ultra-local actions with quantized chiral charges, which are 
therefore doubled (staggered fermions) , and chiral actions which evade the Nielsen- 
Ninomiya theorem, the related cases of domain wall and overlap fermions. 



3.1. Wilson Fermions (undoubled, non-chiral, ultra-local) 

We can alter the dispersion relation so that it has only one low energy solution. The 
other solutions are forced to E ~ l/a and become very heavy as a is taken to zero. 
The simplest version of this solution, called a Wilson fermion, adds an irrelevant 
operator, a second-derivative-like term 

S w = ~ Mipn+n - 2V>„ + VW) - an^V (39) 

to S nalve . The parameter r = 1 is almost always used and is implied when one 
speaks of using "Wilson fermions." The propagator is 

1 — Z7 U sin p u a + ma — rY],, (cos p u a — 1) 

-g(p) = " 2 p \ iry Pm (40) 

a 2^ sm Pn a + i ma ~ r 2^( C0S £V a ~ l )r 

It remains large at ~ (0, 0, 0, 0), but the "doubler modes" are lifted at any fixed 
nonzero r so G(p)^ 1 has one four-component minimum. 

There are actually two dimension-five operators which can be added to a fermion 
action. The Wilson term is just one of them. The other dimension-five term is a 
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Fig. 3. 



The "clover term" 



magnetic moment term 

Ssw ~ -^-i>{x)a tlv F tlv il){x) (41) 

and if both terms are included, their coefficients can be tuned so that there 
are no O(a ) or 0(ag 2 ) lattice artifacts. This action is called the "Sheikholeslami- 
WohlerJ22' or "clover" action because the lattice version of i*} t „ (we'll call it C M „ 
below) is the sum of imaginary parts of the product of links around the paths shown 
in Fig. These days, pure Wilson fermions are rarely used, having been replaced 
by the clover action (with various choices of clover term). 

With Wilson fermions it is conventional to talk about "hopping parameter" 
k — i(ma+4r)~ 1 , and to rescale the fields tp — > v / 2k , 0- The action for an interacting 
theory is conventionally written 

s = E„ ^ni>n - K J2 n ^n(r - 7 M )t/ /1 (n)'0 n+A( + $ n (r + j ^17^-^) 

+c S w^{x)<7^ v C^{x). (42) 

(43) 

From an operational point of view, Wilson-type fermions are closest to the contin- 
uum formulation- there is a four component spinor on every lattice site for every 
color and/or flavor of quark. The simplest bilinears {tp{x)^ ll ^{x), for example), 
could be used as interpolating fields, although in practice more complicated con- 
structs are used to remove lattice artifacts from simulation results. 

Wilson-type fermions contain explicit chiral-symmetry breaking terms. This is 
a source of many bad lattice artifacts. The most obvious is that the zero bare quark 
mass limit is not respected by interactions; the quark mass is additively renormal- 
ized. The value of bare quark mass m q which the pion mass vanishes, is not known 
a priori before beginning a simulation. It must be computed. Simulations with Wil- 
son fermions have to map out vs am q by direct observation. A second serious 
problem with the loss of chiral symmetry breaking is the mixing of operators which 
would not mix in the continuum. This compromises matrix element calculations. 
(See Sec. 6.2 for more discussion.) Finally, nothing prevents the Wilson-Dirac op- 
erator on a gauge configuration from developing a real eigenmode A at any value. If 
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—A happened to equal the bare quark mass dialed into the program, D + m would 
be non-invertible. In practice, these "exceptional configurations" make it difficult (if 
not impossible) to push to small values of the quark mass with non-chiral fermions. 
Rare indeed is the Wilson or clover calculation done at pseudoscalar-to- vector mass 
ratio much below 0.6. 

Besides the simple actions I have described, there are Baroque variations which 
are designed to improve either the kinetic properties (dispersion relation) , or chiral 
properties (minimize additive mass renormalization), or botH^™ 

A recent development!^ which removes exceptional configurations is "twisted 
mass QCD." This is a scheme for Nt = 2 flavors in which the lattice Dirac operator 
is expanded to be 

D twist = D w + ifij 5 r 3 . (44) 

The isospin generator T3 acts in flavor space. The extra term is called the "chirally 
twisted mass." It protects the Dirac operator against exceptional configurations for 
any finite [i: det D tW i S t — det(DyyDw + [J 2 )- It is an amusing exercise to use the 
axial transformations ip — > exp^ar^j^) , ip — ► ip exp^iar^j^) to disentangle what 
would otherwise appear to be an unfortunate mixing of opposite-parity operators. 
I am aware of quenched calculations with this actional but so far, there are no 
published full QCD simulations. 



3.2. Staggered or Kogut-Susskind Fermions (chiral, doubled, 
ultra-local) 

The sixteen-fold degeneracy doublers of naive fermions can be condensed to four 
by the local transformation ip n — > f2„x„, i/> n — > X«OJj where 

3 

There are sixteen different f2's. Using 

fitQ„ = 1; 

fit 7M fi„ +A = (_l)"o+nx+...+» < .- I = «», 

(46) 

we rewrite the action as 

S = y^^»(7p • A m + m)ip n = ^2xn{a(n) ■ A + m)x„. (47) 

n n 

Written in terms of x, the action is diagonal in spinor space. Although we did the 
derivation for free-field theory, it is true for any background gauge field, x is a four- 
component spinor, but since all components interact independently and identically, 
we can reduce the multiplicity of naive fermions by a factor of four simply by 
discarding all but one Dirac component of \. These are "staggered fermions." 
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It is natural to think of the 16 components of a staggered fermion as a fourfold 
replication of four Dirac components. The replication is called "taste." "Taste" is 
the modern word for what used to be called "flavor," as in the sentence "a single 
staggered fermion corresponds to four flavors/tastes." 

There are (at least) two ways to think about the taste content of free staggered 
fermions. The simplest is just to work in momentum space and break up the lattice 
Brillouin zone into 16 components, labeling s M = 0, 1 for each direction, with a 
reduced zone for each component (the "central" one is — n/(2a) < Pfi{s) < n/(2a). 
The massless action is 

S ~ Yl [ ^~P) E ^ sin(2v)V(p) (48) 

The "hypercubic" decomposition is used in simulations. Break the lattice up into 
2 4 site hypercubes and bundle the fields in the hypercube together. Now we define 
the first and second block derivatives (b = 2a) by 

^Xn(N) = ^( X n(N + (i)-Xn(N-fi)) (49) 

D, Xn (N) = ^N + ^) +Xn (N-^)-2 Xn{ N) (5Q) 



In this basis the action is 



X.fl 



b A ip(x) 



ip(x) +mb 4 \^ ip(x)I®Iip(x)(51) 



The sum over x{— Nb) runs over all hypercubes of the blocked lattice. (The notation 
is (spin ® taste). Taste symmetry is four- fold, so the use of Dirac matrices for 
its generators is natural.) The presence of the (75 ® 7^75) term shows that the 
hypercube decomposition still has flavor mixing away from the a — > limit. 

In the continuum, QCD with four degenerate massless flavors has an SU (4)^ (g) 
5i7(4)_R(g) U(l)v symmetry, which spontaneously breaks to SU(4)v, and the pseu- 
doscalar spectrum consists of 15 Goldstone particles and a massive (from the 
anomaly) eta-prime. On the lattice, taste and spin rotations are replaced by shifts 
and rotations in the hypercube. The continuous symmetries of the continuum are 
broken down to discrete symmetries. In particular, taste symmetry is broken. Only 
a U{l)v <8> U(1)a survives, and the spontaneous breaking of the lattice analog of the 
flavor non-singlet axial symmetry produces a single Goldstone boson at nonzero a. 
The other "would-be" Goldstones are non-degenerate pseudoscalar states whose 
mass goes to zero in the continuum limit. (QCD with Nf flavors of staggered 
fermions has an internal SU(Nf)i J ® SU(Nf)n chiral symmetry.) The U(l) chiral 
symmetry protects the quark mass from additive renormalization, and staggered 
fermions are preferred over Wilson ones in situations in which the chiral properties 
of the fermions dominate the dynamics. 

There are two complementary ways to think about taste breaking. In momentum 
space, taste mixing occurs in perturbation theory via the emission and absorption of 
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Fig. 4. An example of flavor symmetry breaking in an improved staggered action. The different 
7's are a code for the various pseudoscalar states. Data are from Ref. 1151 For an explanation of 
the splitting, see Ref. 1261 



gluons with momentum near — tt/cl. These kick quarks from one s M momentum 
sector into another one. In the hypercube basis, different tastes sit on different 
locations in the hypercube. The environment of link variables varies across the 
hypercube, so different tastes see different local gauge fields. 

The pseudoscalar spectrum shows an interesting (approximate) degeneracy (see 
Fig. 21 for an example), first described by Sharpe and Led^l It is reminiscent of the 
pattern of splittings of the energy levels of an atom in a crystal lattice. 

Today's approach to staggered fermion simulations exploits the conversion of a 
U(l) chiral symmetry for a single staggered flavor (with its single would be Gold- 
stone boson) into the SU(Nf) chiral symmetry of Nf flavors (with N* — 1 Gold- 
stones). QCD with three flavors of quarks is treated with a lattice model with three 
flavors of staggered fermions, each with its own mass. Each flavor would correspond 
to four tastes. Three of the four tastes (per flavor) have to be excised from one's 
predictions. How that is done is different for valence and sea quarks, and will be 
described below. 

This decomposition has been around since the mid-90 'J2ZESI 

It was introduced 

in order to understand weak matrix elements of staggered fermions. Recent tech- 
niques for staggered fermions, which make heavy use of chiral Lagrangians, have 
brought it to greater importance. 

Taste symmetry violation for a single staggered flavor basically means that 
interactions mix the different tastes. For the connection between staggered flavor 
and physical flavor to work, taste symmetry violation must be minimized. People 
attack this problem in two ways. 
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First, they modify the action to suppress taste violations. This is commonly done 
by replacing the thin links of the fermion connection by a fat link. The "Asqtad" 
fat link actional (which also adds a third - nearest - neighbor coupling to improve 
the dispersion relation) has seen the most extensive use. Other choice^^l reduce 
taste violations more. 

Second, lattice data is analyzed including the effects of taste violations. The key 
to doing this was provided by the Sharpe and LeJ^ analysis of taste mixing, whose 
construction of a low energy chiral effective theory including explicit taste-breaking 
interactions predicted t he degeneracies shown in Fig. 01 Their work was generalized 
by Aubin and Bernard^ to Nf flavors of staggered fermions (4 Nf tastes). One 
introduces a field £ = exp(i$//), a ANj x 4 Nf matrix, and $ is given by: 



4> 



/ U 7T+ K+ 
7T" D K° 

K- K° S 

V ; ; 



/ 



(52) 



The entries are the meson fields composed of different staggered flavors, and are 
traces over the 16 taste-product qq bilinears of each staggered flavor. The mass 
matrix is 



M 



(m u I •• 
m d I •• 
m s I ■ ■ 

V i ••/ 



(53) 



The Lagrangian is 



o 4 
+ 2 -^ {Ul + D I + S I + ---) 2 +a 2 V, 



(54) 
(55) 



where the term weighs the analog of the flavor singlet rf . (The "J" subscripts 
display that this involves the taste singlet term for each flavor.) The a 2 V term is 
the taste-breaking interaction, a sum of terms quadratic in E with various taste 
projectors, parameterized by six coefficients (only one is big). Now one computes 
"any" desired quantity with this Lagrangian, typically to one loop, as a function 
of quark masses and all other coefficients. Parameters of Nature are determined 
when mass-dependent Monte Carlo data is fit to this functional form. For example, 
a one-loop fit to mps/(fni + 7712) for the true would-be Goldstone boson made of 
quarks of mass mi and ni2 would involve /x, /, two Gasser-Leutwyler parameters, 
three otherwise unconstrained lattice parameters, and involves chiral logarithms 
whose arguments are all the observed pseudoscalar masses. Fits to or fjc are 
similar. 
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It is crucial to analyze staggered data this way. In 2001 the MILC collaboration 
tried and and failed to fit their pseudoscalar masses to the continuum chiral loga- 
rithm formula. The same data, with fits which include taste violations, forms part 
of their recent high precision calculation of hadronic parameters (see Sec. 5). 

3.3. Evading the no-go theorem - domain wall and overlap 
fermions 

Two related schemes allow one t o evad e the Nielsen-Ninomiya theorem. The first 
action is the "domain wall fermior PM^ !" It is a variation on the idea that a fermion 
coupled to a scalar field which interpolates between two minima (a soliton) will 
develop a zero-energy chiral mode bound to the center of the soliton. Embed QCD 
is a five-dimensional brane world with a kink. This is a discretized fifth dimension. 
The gauge fields remain four dimensional. A chiral fermion will sit on the four 
dimensional face of the kink. This is the chiral "domain wall fermion." As Kaplan 
puts it, "The extra dimension is the loophole in the Nielsen-Ninomiya theorem 
through which the fermions have wriggled." 

Several groups have well-developed programs of QCD phenomenology, both 
quenched and dynamical, with domain wall fermiona^. We will encounter their 
results in the section of this review dealing with kaon physics. 

The practical consideration which must be dealt with in domain wall simulations 
is that in the computer, the fifth dimension is not infinite. There will be an anti- 
kink somewhere else in the fifth dimension or some equivalent boundary surface, 
with an opposite-chirality fermion pinned to it. As long as the two kinks are far 
away, the fermionic chiral mode on the kink doesn't see the mode on the anti-kink 
and the 4-d theory on the kink will remain chiral. But if the anti-kink is too close 
(typically, if the fifth dimension is too small or the fermion eigenmodes on the 
kink are insufficiently localized in the fifth dimension) the modes mix and chiral 
symmetry is broken. The fermion will be observed to pick up a small additive mass 
renormalization. How close is "too close" is (yet) another engineering question, and 
is dealt with in the usual way (by modifying lattice discretizations). 

The four-dimensional analog of this formulation uses the Ginsparg- Wilson rela- 



Its explicit realization is Neuberger's "overlap fermion- -<" The massless overlap 
Dirac operator is 




j 5 D(0) + £(0) 75 



—D(0) 15 D(0). 



(56) 



D(0) =x (l + 



z 



) 



(57) 




where z = d(—Xo)/xo = (d— xq)/xq and d(m) = d + m is a massive "kernel" Dirac 
operator for mass m. d can be (almost) any undoubted lattice Dirac operator. The 
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chiral symmetry of these remarkable actions prevents additive mass renormaliza- 
tion, preserves all continuum Ward identities (up to contact terms) and knows 
about the index theoreirP^. 

The eigenmodes of D(0) sit on a circle of radius xq centered at (xq,0) in the 
complex plane. The zero mode is chiral (as is the mode at (2xq, 0)) and the complex 
modes are non-chiral and paired (with complex conjugate eigenvalues). The Nielsen- 
Ninimiya theorem is evaded because the chiral charges are not quantized. One 
waySDto gee ^jjjg j s ^ rea ij ze that the action S = -0£>(O)-0 is invariant under the 
gauge field dependent axial transformatiorP^ 

Sip = T%ip = Tj 5 (l - —)tp\ Sip = -075T (58) 
x 

where T is a U(Nf) flavor rotation generator. Directly from the Ginsparg- Wilson 
relation, (75) 2 = 1, so the combination of this axial vector transformation, plus the 
usual vector one, generates a conventional current algebra. The fermion measure 
is not invariant under the U(1)a version of Eq. EH1 and this leads to the anomaly. 
For small D, Eq. EHlis the usual chiral rotation. This is the situation for the low 
eigenmodes of the Dirac operator. But at the other corners of the Brillouin zone 
D ~ 2xq, 75 flips sign, and the transformation does not correspond to a chiral 
rotation. One can also consider a symmetric version of Eq. 1581 

# = T 75 (1-^)V; ^ = ^(1-^)T 75 (59) 

At the "far corners" of the Brillouin zone, where D = 2xq, this transformation 
vanishes. Either way, the eigenvalue of the chiral charge is not quantized, it varies 
with D. 

The massive overlap Dirac operator is conventionally defined to be 

D(m q ) = (l-^-)D(0) + m q (60) 

and it is also conventional to define the propagator so that the chiral modes at 
A = 2xo are projected out, 

D-\m q ) = - l ——{D-\m q )~^-). (61) 

l-m, q /(2x ) 2xq 



This also converts local currents i nto order a 2 improved operator^ 

No overlap action is ultra-loca^S . Satisfying the Ginsparg- Wilson relation re- 
quires an action which has connections spread out to an arbitrary number of lat- 
tice spacings. For sufficiently smooth gauge configur atio ns, Hernandez, Jansen, and 
Luscher have proved that overlap actions are locai^l Golterman and ShamiiEl 
have presented a convincing argument that sufficiently rough gauge configurations 
can drive overlap and domain wall fermion actions nonlocal. It is presently an item 
of debate, how much today's simulations are corrupted by non-locality. 

The hard part of an overlap calculation is the "step function" (e(z) = j 5 z/Vz^z). 
There are various tricks for evaluating it, basically as polynomials in z (Chebyshev 



1, 2008 22:50 WSPC/INSTRUCTION FILE draftl 



Lattice QCD 23 

polynomials) or as a ratio of polynomials A{z) / B{z) ~ r )%z'^ l {\/{z*z 4- c n ). (For 
discussions of overlap technology, see Refs. 03031) The degree of success of these 
evaluations typically depends on the conditioning number of z, and to improve 
that, people remove low eigenmodes of z from the evaluation and treat them ex- 
actly. Overlap calculations typically cost a factor of 50-100 as much as ordinary 
Wilson calculations, but to quote that fact alone is unfair: they can be used to 
study small quark mass (chiral) physics with full flavor symmetry at quark masses 
where Wilson-type actions simply fail due to exceptional configurations. Unfortu- 
nately, overlap fermions are still too expensive for anything but extremely tentative 
dynamical fermion simulations. 

There is of course an aesthetic (or engineering) debate among lattice practi- 
tioners over whether it is better to do simulations which have approximate chiral 
symmetry or exact chiral symmetry. The former are generally computationally less 
expensive, but they have to be checked against the appearance of unwanted effects. 
With an exact algorithm, one might sleep better at night knowing that the calcu- 
lation does not have chiral artifacts, but it might also happen that one is simply 
unable to generate an interesting data set for analysis. 

Because simulations with domain wall and overlap fermions are so expensive, 
there is a tendency to cut corners, either by making the lattice spacing uncomfort- 
ably large, or the simulation volume too small. (The author's calculation of Bk 
using overlap fermions is done at a small enough volume, that baryon masses are 
clearly affected-they are pushed to artificially high mass by the squeezing of the 
box.) 

One can turn a bug into a feature by simulating QCD in the "epsilon regimd 
This is QCD in a box whose length L is small compared to the pion Compton 
wavelength, m^L << 1. It is not small simulation volume, since L should also be 
large compared to any 'typical" QCD confinement scale. Thi s ye ar we have begun 
to see simulations (mostly with overlap fermions) in this regime^! By analyzing the 
behavior of hadron correlators, and matching onto chiral Lagrangian calculations, 
also done in the epsilon regime, low energy properties of QCD may be extracted. 

The descriptions of domain wall and overlap fermions did not make them look 
particularly similar. The connection comes when one takes domain wall fermions 
and integrates out the bu lk field s to construct a four dimensional effective action 
of the light "pinned" fieldJ™, 

(1 + H) Nb - (1 - H) Ns 

where H — 752/(2 + 2:) and z was defined above. As the length of the fifth dimension 
becomes infinite and the lattice spacing in the fifth dimension vanishes, this Dirac 
operator becomes the overlap operator. 

Finally, different kinds of fermions have different constraints on the number of 
flavors which can be simulated with various algorithms. Consider a single flavor of 
staggered (four tastes) or Wilson-type fermions (one flavor), with a massive Dirac 
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operator M. In either case M' — 75M75. The functional integral for a theory 
with a single flavor is just det M. The massless staggered fermion Dirac operator 
M{m = 0) is antihermetian so its eigenvalues lie along the imaginary axis, and they 
come in pairs For massive fermions, det M is a product of (iX + m)(—i\ + m) 
factors, thus it is positive-definite for all gauge configurations. Standard algorithms 
(Refreshed Molecular Dynamics®, Hybrid Monte Carlo®) actually work with M^M, 
which would redouble the number of flavors. But because M^M only connects 
"even" sites on the lattice with other "even" sites (and 'odd" with "odd"), one 
can take the square root of the operator by removing the fermion fields on all the 
"even" (or "odd") sites. 

Wilson fermions have eigenvalues which are real, and complex conjugate pairs. 
The determinant from the complex conjugate pairs is basically just like the stag- 
gered determinant; it is a product of (m + A r ) 2 + A 2 . This part is positive definite. 
However, nothing protects the real eigenvalues from taking any value, so the prod- 
uct of real eigenvalues could have any sign, and could change sign as the gauge 
field evolves. With two degenerate flavors (or Nf — 2j), one can get a positive 
determinant because detM = detM 1, and (detM) 2 = det(M t M). 

Furman and Shamir have proved that one can simulate any number of flavors 
with domain wall fermions. For overlap fermions, the real eigenvalues are restricted 
to be cither zero or a positive number, and so one could also simulate any number 
of fermions, each one with its own detM, with no problem with positivitjl^l 



4. The Quenched Approximation in Health and Sickness 

For many years, calculations of spectroscopy and matrix elements in QCD have 
used the "quenched approximation" in which the fermion determinant is completely 
discarded from the functional integral. This approximation was adopted pretty 
much for reasons of expediency, because numerical simulations in full QCD were 
simply too expensive. From its earliest days it has produced results for known 
quantities which were in remark able agreement with experiment. The results of one 
very careful quenched simulation^! f the light hadron spectrum are shown in Fig. 
El One could not ask for anything better. 

The quenched approximation has many of the ingredients of successful hadron 
phenomenology. Quarks are confined (with a linear confining potential if they are 
heavy). Chiral symmetry is spontaneously broken. In it, all states are (at first 
glance) infinitely narrow, because qq pairs cannot pop out of the vacuum. One 
might also try to "justify" the quenched approximation by an appeal to the quark 
model: in the quenched approximation, all mesons are qq pairs, and all baryons are 
qqq states. This also appears to be rather similar to the large- N c limit of QCD. 

However, the situation has changed. For about ten years people have identified 
specific situations where the quenched approximation would give qualitatively dif- 
ferent behavior than QCD with any nonzero number of sea quarks. This behavior 
now needs to be dealt with - or is beginning to be seen - in simulations. 
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Fig. 5. Comparison of quenched results from Ref, 1481 with experiment. 



The best way to see what is going on is to consider the low energy limit of 
QCD, and not to think about quarks and gluons, but in terms of an effective 
field theory theory of QCD, described by chiral Lagrangians in which the would- 
be Goldstone bosons are fundamental fields. These Lagrangians have a set of bare 
parameters (quark masses, f v , the quark condensate S, . . .). As far as the chiral 
Lagrangian is concerned, these are fundamental parameters. As far as QCD is 
concerned, one could compute these parameters from first principles (for example, 
/^TOtt = (O|V'7o75'0l 7r ))i this would fix the parameters of the chiral Lagrangian, 
and then one could throw away the lattice and compute low energy physics using 
the chiral Lagrangian. Quenched QCD and QCD with nonzero flavor numbers are 
different theories and their low energy parameters will be different. But there is 
more. In full QCD the eta prime is heavy and can be decoupled from the interactions 
of the ordinary Goldstone bosons. In quenched QCD the eta prime is not really a 
particle. The would-be eta prime gives rise to "hairpin insertions" which pollute 
essentially all predictions. 

Let's consider the eta prime channel in full QCD and quenched QCD in various 
extreme limits. In ordinary QCD, the eta prime propagator includes a series of 
terms in which the flavor singlet qq pair annihilates into some quarkless state, then 
reappears, over and over. This is shown in Fig. EI Let's assume we have N v valence 
quarks, N s sea quarks, and N c colors. (The case N v ^ N s with N s ^ is called 
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Fig. 6. The eta-prime propagator in terms of a set of annihilation graphs summing into a 
geometric series to shift the eta-prime mass away from the mass of the flavor non-singlet pseu- 
doscalar mesons. In the quenched approximation, only the first two terms in the series survive as 
the "direct" and "hairpin" graphs. 

"partial quenching" as opposed to the N s — quenched approximation.) Let's sum 
the geometric series for the eta prime propagator 

V '(q) = C(q)-H (q)+H 1 (q) + ... (63) 

where C(q) = 1/d, d = q 2 + m 2 , is the "connected" meson propagator, the same 
as for any other Goldstone boson. H n is the nth order hairpin (with n internal 
fermion loops). Simple N c counting appropriate to the large- N c limit (taking g 2 N c 
fixed), with flavor and color singlet sources and sinks, shows that C(q) oc N®Ny 
and H n (q) oc 1/N?. 

The lowest order hairpin is 

= < M > 

where A 2 is 0(N°). Each extra loop gives another factor of —N S X 2 / (N c d). Summing 
the geometric series, we find 

^ q) = d {1 -wJ + W s dTx^jN c - (65) 

This formula has a number of interesting limits. First, if N v = N s = Nf ^ 0, 
we see that we have generated a massive rj propagator with m 2 , = X 2 Nf/N c . 

The next case is N v < N s . There we have both the Goldstone mode and the 
eta prime propagator. This is also as expected: think of the SU (3) flavor symmetry 
limit and imagine a source uu + dd: it couples to a mixture of the eta (a Goldstone) 
and the eta prime. 

Now we could take iV c to infinity at fixed Nf. With iV c -scaling of the vertex, 
the eta-prime mass falls to zero, and it becomes a ninth Goldstone boson as U(1)a 
is restored. At any finite N c , the eta prime is still an ordinary particle. 

However, the quenched limit is different-it is N s = first. In that case (of 
course) 

"® = d-dN- e X d' (66) 
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In the eta prime channel there is an ordinary (but flavor singlet) Goldstone boson 
and a new contribution-a double-pole ghost (negative norm) state. In the N c — oo 
limit, the double pole decouples, but finite N c quenched QCD remains different 
from finite- N c full QCD. The limits of large N c and quenching don't commute. 

Similar behavior will occur in any other flavor singlet channel (think about the 
ijj meson), of course. Lattice people haven't talked about it because the signals are 
noisier, and because there is no chiral Lagrangian paradigm. 

The double pole would reappear if the sea quarks and valence quarks had dif- 
ferent masses. Its residue would be proportional to the difference of the sea and 
valence pseudoscalar squared masses. A convenient form of partial quenching is in 
fact to compute hadronic properties for one set of sea quark masses (because each 
set is expensive) and many values of valence quark masses (because each set is 
cheap), 

(What is the internal flavor symmetry group of partially quenched QCD? The 
determinant of the valence quarks is not present in the functional integral. One 
can cancel the determinant by introducing N v valence bosonic quarks, s ince th e 
power of the determinant is negative for bosons. The flavor group become^^Ola 
graded group SU(N V + N s \N v )l ® SU(N V + N s \N v )p which spontaneously breaks 
to SU(N V + N S \N V ). The low energy effective Lagrangian has additional mesons, 
corresponding to bound states of quarks and bosonic quarks.) 

Whe re the eta-prime comes in is in the calculation of corrections to tree-level 
relationJ23El These are typically dominated by processes with internal Goldstone 
boson loops, contributing terms like 



(plus cutoff effects). The eta-prime hairpin can appear in these loops, replacing 
G(£:, to) — > —G(k, m)X G(k, to) and altering the chiral logarithm. Thus, in a typical 
observable, with a small mass expansion 



quenched and Nf — 3 QCD can have different coefficients (different -B's in Ea. 168(1 . 
seemingly randomly different. (Quenched /„■ has no chiral logarithm while it does 
in full QCD, the coefficients of 0+, the operator measured for Bk, are identical in 
quenched and full QCD, etc.) Even worse, one can find a different functional form. 
For example, the relation between pseudoscalar mass and quark mass in full QCD 
is 




(67) 



Q(tops) = A(l+B 




logmps) + . . . 



(68) 



mpg = Am q (l + 




Llog(m 2 /A 2 )] + .... 



(69) 



In quenched QCD, the analogous relation is 



(mps)7K) = A[l - <5(ln(TO 2 /A 2 ) + 1)] + 



(70) 
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where S — \ 2 /(8ir 2 N c f 2 ) is expected to be about 0.2 using the physical rf mass. 
This means that rn 2 PS /m q actually diverges in the chiral limit! Many quenched 
simulations actually search for these "quenched chiral logarithms," and a 
claim to have seen them. 

Quenched QCD will not go away any time soon, but its days as a source of 
precision numbers for QCD matrix elements are clearly numbered. It will continue 
to be used for tests of methodology as a proving ground for new ways of processing 
data, and a way to settle some controversies. There is also quite a bit of contin- 
uum phenomenology which can be altered to appl y to quenched QCD and which 
simulations of quenched QCD can validate - or noll^l 

One example of a controversy quenched QCD can address is the following: At 
Lattice 2000, S. AokP2l argued that the continuum limit of spectroscopy with stag- 
gered quarks and Wilson- type quarks might be different. This would represent a 
loss of universality, and would represent a serious problem for (at least) one kind 
of lattice discretization of fermion. The data he showed used unimproved Wilson 
fermions (which have order a lattice artifacts), unimproved staggered quarks, and 
clover quarks (both of which have 0{a 2 ) scaling violations). His observed disagree- 
ment may have been because of the large lattice artifacts in (at least some) of the 
data sets, which make extrapolations to a — > difficult. An example of more recent 
data is shown in Fig. This is an example of an "APE plot" which compares two 
dimensionless ratios (in this case (nips /my) 2 vs /my)- If scaling violations are 
absent, data from different lattice spacings will lie on a common curve. The data 
shown in Fig. are not extrapolated to the continuum limit, but they do give rea- 
sonable evidence that there is not a lot of scale violation in quenched spectroscopy 
done with improved actions. 



5. Simulations with Dynamical Fermions 

Simulations which include dynamical fermions go back about fifteen years. Because 
they are so costly, people begin with fairly heavy sea quarks (pseudoscalar/vector 
meson mass ratios > 0.7), and try to go down. Volumes tend to be smaller, and lat- 
tice spacings tend to be coarser than quenched simulations. Statistics are a problem; 
the simulations have long autocorrelation times. 

Over the years I can recall simulations with 2, 3, and 4 flavors of fermions, using 
staggered, Wilson, and clover quarks. The RBC collaboration is in the earliest stage 
of a dynamical fermion simulation of domain wall fermionJS3. The first simulations 
with overlap fermions have just appeared: they cost easily a factor of a hundred 
times ones with Wilson quarks^l The lattice conferences have had a steady diet 
of spectroscopy and matrix element calculations fr om these projects. But this year 
there seemed to be a pause. I think many people^l have decided that they just 
cannot push the quark masses down far enough to be interesting, and have gone 
back to studying algor (I expect to get a lot of unhappy mail about this 

sentence.) 
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Fig. 7. A comparison of staggered spectroscopy, using an improved action, (pluses and octagons) 
at two lattice spacings (from Ref . 55]) with that from an improved Wilson- type action (from Ref . 
1561 in crosses. 



There is one notable exception: The most interesting lattice calculations with 
dynamical simulations are the ones being done by the MILC collaboration. They 
have combined with other groups to do calculations of spectroscopy with three 
flavors of dynamical staggered quarks, a strange quark at approximately its cor- 
rect physical value, and degenerate up and down quarks whose masses run down 
to m s /5. Their gauge configurations are used as backgrounds for simulations with 
heavy flavors. The analysis uses the taste-breaking chiral Lagrangian described in 
Sec. 3. They have data at two lattice spacings, 0.13 and 0.09 fm, large volumes, 
and high statistics. They have studied a variety of processes for which lattice meth- 
ods ought to work well: spectroscopy of hadronic states which are not close to 
decay thresholds, properties of pseudoscalar mesons, hadrons with strange valence 
quarks. Their recent preprinlP^ presents a remarkable agreement of simulation with 
experiment, shown in Fig. El 

The MILC part of the staggered project is several years old. It has produced a 
lot of interesting spectroscopy. One can clearly see the effects of sea quarks in the 
data. Here are some exa mplefPl 

The heavy quark potential in QCD looks something like V(r) ~ err — C jr. 
The location of the bend in the potential is related to a quantity called the Sommcr 
parameter ro (see Ref. l63[) via the force, —r 2 F(r) = 1.65 at r — r$. This corresponds 
to ro of about 0.5 fm. Fig.|5|shows the dimensionless quantity r^^fa as a function of 
the quark mass, represented by (m w /m p ) 2 . This places the quenched approximation 
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Fig. 8. Comparison of quenched results with results from simulations with 2+1 flavors of stag 
gered fermions, from Ref. 1611 



at (to^/top) 2 = 1, and the chiral limit at the left side of the graph. In these plots the 
octagons are runs with three degenerate sea quarks, except for the rightmost point 
which is the quenched limit. Squares are runs with am s = 0.05 at a = 0.13 fm. This 
is roughly the physical value of the strange quark mass. For these runs, am Ui d < 
0.05. The isolated diamond is a two flavor run. Finally, the cross at (m„-/m p ) 2 = 1 is 
an 0.09 fm lattice spacing quenched run. From the two quenched points the authors 
infer that remaining lattice artifacts are small compared with the effects of the sea 
quarks. One can clearly see the distinction between two and three flavors, as well 
as the effect of using two light and one heavy flavor rather than three degenerate 
flavors (the "kink' at (rnjrn p ) 2 w 0.46). 

Lacock and Michael*^ have observed differences between the quenched meson 
spectrum and the real world. They studied the quantity 

dray 



J = rriK* 



dm 2 PS 



(71) 



where my and mps are the vector and pseudoscalar meson masses. This quantity 
has the advantage of being relatively insensitive to the quark masses, so that accu- 
rate tuning of the strange quark mass or extrapolation of the masses to the chiral 
limit is not essential. 

Of course, to compare to experiment the derivative in this expression must be 
replaced by a ratio of mass differences, and MILC choose 



J = rriK* 



2{mj 



(72) 



Here m p is the mass of the vector meson including two light quarks, etc. Figure ITUI 
shows the results for J in quenched and three flavor QCD. This is plotted versus 
vtlk* /rriK, for which the real world value is 1.8. The burst is the real world value 
of this definition of J (0.49), and the cross is the value of J found in the quenched 
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Fig. 9. Effects of dynamical quarks on the shape of the potential: ro^/a as a function of the 
quark mass. The two quenched points are at the far right, with the octagon coming from the 0.13 
fm run and the cross from the 0.09 fm run. The remaining octagons are full QCD runs with three 
degenerate flavors, and the squares are full QCD runs with two light flavors and one heavy. The 
diamond is the two flavor run, and the burst at the left is a linear extrapolation of the 2 + 1 results 
to the physical value of (m„- /m p ) 2 . From Ref. 1621 

simulations of Ref. E] One can see a clear effect of the sea quarks on this quantity. 
Figure HTH also contains one point with two dynamical flavors. This point falls near 
the three flavor points, indicating that the dynamical strange quark is less important 
than the two light quarks. 

If there are dynamical fermions in the configurations, then one ought to see pro- 
cesses like particle decay in the simulations. What the lattice measures is energies. 
If a particle A could decay into a BC state, the A — A correlator would have a 
contribution which would show an exp(— Esct) behavior, with Esc — m B + rnc- 
Ea and Esc might depend on simulation parameters (like quark mas ses) and one 
might look for avoided level crossings as the parameters were varie (Spin mod- 
els with several species of excitation show this kind of behavioral) The chief target 
is the rho meson. Evidence for a pint coupling has not been seen. The problem is 
that the rho is J = 1, so in the decay p — > 7T7t the pions are P-waves. The mini- 
mum energy of a two-pion state is then + y/ m\ + pf nin , where p m in = 2tt/L is 
the smallest nonzero momentum in a box of width L. This effectively pushes the 
threshold for the level crossing to smaller quark mass. (For further discussion of 
quantum number effects for staggered fermions, see Ref. 67 ) 

However, the instability of an excited state has been seen in a different channel, 
the isotriplet scalar clq mesor^. It is clearly very different in the quenched and full 
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Fig. 10. The mass ratio "J" in the quenched and full QCD calculations. Squares are the three 
flavor results, and octagons are the quenched results. The diamond is the two flavor run, using a 
non-dynamical strange quark with mass am q = 0.05. The burst is the real world value, and the 
cross is the quenched value of Ref. 64 The smaller error bar on the cross is the statistical error, 
and the larger the quoted systematic error. From Ref. 1551 



QCD runs. For large quark masses there is no visible difference, but as the quark 
mass is decreased the full QCD ++ mass drops below all the other masses. For all 
but the lowest quark mass, the quenched ++ is close to the other P-wave meson 
masses. It is plausible to ascribe the behavior of the full QCD mass to the decay 
of the eto into ir + r\. (Bose symmetry plus isospin forbids decay into two pions.) 
Figure ITU illustrates this interpretation. In the figure I plot the quenched and full 
++ masses versus quark mass. The straight line in the graph is a fit to the quenched 
mass for the heavier quarks, and represents the mass of a qq state. The curved line 
with the kink at am q = 0.05 represents the mass of it + r\. For am q > 0.05 MILC 
used three degenerate quark flavors, so the r\ and 7r are degenerate and this line is 
simply twice the pion mass. For am q < 0.05 one does not have direct information 
on the r\ mass, so the Gell-Man-Okubo formula written in terms of an "unmixed 
ss" mass (just the pseudoscalar mass at am q = 0.05) is used: 

m\ = {ml + 2tt4)/3. (73) 

In the quenched case the ao mesons can couple to two-meson states through a 
"hairpin diagram" on one of the meson lines. Such diagrams can behave like powers 
of t times e -2 ™' 1 and therefore masquerade as a light ao when 2m 7r < m ao . This 
may explain the lightest quark mass quenched poiirJEsl. 

The authors of Ref. HTJJ argue that their results show that at last, there is a 
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Fig. 11. + + masses versus quark mass. The lightest fitted energy in the scalar channel. Octagons 
are quenched results, squares are three flavor results, and the burst is the two flavor run. The 
straight line is a crude extrapolation of the heavy quark points. The curved line is the it + 77 mass 
estimate, as discussed in the text. The short vertical line marks the approximate quark mass where 
the ao mass is twice the quenched pion mass. From Ref. 62 

workable algorithm for simulating QCD with light fermions, and that its successes 
mean that the time is ripe to apply it to a wide variety of lattice matrix element 
calculations. 

However, this conclusion is not universally accepted by the lattice community. 
Many people are worried about the way staggered fermions are used to represent 
2+1 flavors of sea quarks. 

Lattice simulations treat valence quarks and sea quarks differently. Basically, 
valence quarks (the ones attached to the external sources) do not have any quan- 
tum numbers (other than their mass and spin). One computes classes of Feynman 
diagrams on the lattice, then re-weights them with global symmetry indices and 
bundles them together. (For example, the same propagators are used for the quark 
and the antiquark in a mass-degenerate meson). In staggered fermions, one uses a 
single flavor of staggered fermions, with its four tastes, and computes correlation 
functions in which the sources project (nearly) onto the same initial and final taste. 
The quark could hop temporarily into a different taste state as it propagates across 
the lattice (this would happen by emitting and absorbing hard gluons), but this is 
just cutoff scale physics which contributes 0(a 2 g 2 ) scale violations. 

But sea quarks are different - the number of flavors matters. In the absence of 
flavor symmetry breaking, mass-degenerate fermions have identical spectra and one 
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can simulate Nf flavors of degenerate fermions simply by raising det M for a single 
flavor to the Nf power. But each taste of staggered fermions contributes to the 
functional integral like a single flavor of continuum quark. Det M is the weighting 
for four continuum flavors. To get the weighting for a single flavor, people doing 
simulations with staggered fermions re-weight the determinant to det 1 / 4 M per 
staggered flavor. 

Are simulations with staggered fermions with det 1 / 4 M fundamentally correct? 
Until recently, all the theory of the fractional d eter minant was a single sentence 
in a 1981 paper by Marinari, Parisi, and RebbP^l "On the lattice, this action 
(Sg — (l/4)Tr InM) will produce a violation of fundamental axioms, but we expect 
the violation to disappear in the continuum limit and then recover the theory with 
a simple fermion." In the past year, the situation has heated up, although the 
literature is still spars probably because there are no well-defined questions. 

It is a peculiar exchange-of-limits question. Any N-flavor action action with fla- 
vor symmetry can be reduced to a one-flavor action by taking the 1/iVth. root. But 
in staggered fermions taste symmetry is broken by gluon interactions at nonzero lat- 
tice spacing. What does taking the fractional root before restoring flavor symmetry 
do? 

Here is my own condensation of arguments I have participated i 

The most serious problem could be a loss of locality. Is there a single-flavor 
(undoubled) local fermion action whose determinant is equal to det 1 / 4 M (or equal 
to it up to cutoff effects)? If such an action exists, then there are no problems 
in principle with det 1 / 4 M staggered fermions. This would certainly be something 
novel! Its chiral symmetry properties would certainly be very different from those 
of ordinary staggered fermions. Its chiral charge could not be discrete, because that 
would violate the Nielsen-Ninomiya theorem. There is no reason that the action 
could not be ultra-local, but it seems hard to combine an ultra-local action with 
an unquantized chiral charge. It's easier for me to imagine that it would have 
couplings falling exponentially with distance. If such an action does not exist, then 
universality is lost. It would not possible to integrate out short distance physics to 
leave a long distance effective action whose couplings flow to a fixed point, because 
remnants of the original action will be left behind after each blocking step. 

And of course, there are many ways to define matrices whose determinants are 
equal. Recall the discussion at the end of the last section about taste constraints. 
Det M was converted into det(M^M) 1 / 2 by deleting fields on half the sites of the 
lattice. That gives a "square root" action which is ultra-local, as opposed to working 
with an action (Ai^M) 1 / 2 , which probably isn't locals 

In the momentum-space decomposition of tastes (Eq. I48JI . each taste occupies 
a segment of the Brillouin zone equivalent to — ir/2 < < tt/2 with a free field 
(27^ sin (p^) dispersion relation. The Fourier transform of its inverse propagator is 
an action with long range (power law falloff) couplings. In the hypercube decom- 
position, each taste's action (the first term of Ea. 151(1 is ultra-local, but tastes still 
mix for any finite value of a through the second term of Eq. 1511 
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It is certainly the case that an action with det 1 ^ 4 M fermions is not unitarity, 
in the sense that absorptive parts of processes are not related to cross sections. 
This is a generic feature of any partially-quenched theory, where the actions of 
valence and sea quarks are different. Different actions mean different Feynman rules. 
(Even though there may be no action corresponding to det 1//4 M fermions, one can 
expand (l/4)TrlnM in a power series in g and construct graphs.) Unitarity relates 
processes in which quarks are internal lines (with internal vertices, propagators, 
etc.) to other processes where the same quarks appear as external lines. Then 
they would be valence quarks, with valence Feynman rules. For staggered fermions, 
the differences occur in processes where taste is changed. There are 4 identical 
vertices (per staggered flavor) which do not change taste in the determinant, but 
only one process which changes taste i to another taste j — 1 . . .4). But taste 
changing interactions involve hard gluons, with momentum transfer ~ tt/&, so 
maybe we are just talking about cutoff effects, again. It is also unknown whether 
there are practical consequences for Euclidean-space lattice calculations of non- 
unitarity: unitarity describes what happens when scattering amplitudes are rotated 
from Euclidean space to Minkowski space. 

It is generally believed that "partially quenched" QCD is well-enough behaved 
to make predictions. Indeed there is an extensive literature^! on the use of partially 
quenched chiral perturbation theory to measure the parameters of full QCD. If the 
det 1/ /4 _M theory is local, then it is just some funny kind of partially quenched 

qceP 

and again nothing would be wrong in principle. 
There are several directions this controversy could play out: 
First, someone could show analytically that there is some local action which has 
the same determinant as det 1 / 4 M, or prove that no such action exists. If such an 
action exists, then there is nothing wrong in principle with det 1 ^ 4 M weighting. It 
would be an engineering question to ask whether this formalism is good enough to 
produce reliable continuum numbers. If no action exists, then QCD with det 1 ^ 4 M 
weighting is not in the same universality class as continuum QCD. It would still be 
the case that its predictions are quite close to experiment (I am not disputing the 
results of Ref. I61f) for many quantities, but one might be able to identify processes 
where it would fail. It would be like the quenched approximation all over again, 
only better: fundamentally uncontrolled, phenomenologically successful. I suppose 
an intermediate possibility is that for sufficiently smooth background configurations, 
det 1 / 4 M could correspond to a local action, but for rough fields it might not. This 
is easy to visualize when one thinks of the fractional power as basically counting 
the geometric mean of the different taste eigenmodes: (Y[f Aj) 1/W . At sufficiently 
rough coupling, where the chiral symmetry is only a U(l), the eigenmodes are no 
longer even approximately degenerate, and the geometric mean would not have 
much meaning. 

Second, one could continue to compute known quantities, see that the simula- 
tions continued to agree with experiment, and go on to do make yet more predic- 
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tions. This is how comparison of theory and experiment generally happens. The 
dangerous possibility here is that something might not agree with experiment, and 
it might be a matrix element from a Standard Model test. Then, is the Standard 
Model wrong, or is it the simulation with a fractional power? 

Checks with simpler models would be very useful. The Schwinger model (for 
different flavor number) is a good test case, and it is simple enough th at it can 
be overwhelmed by numerical simulations. Recently, Diirr and HoelblingZ^I have 
begun a study of the Schwinger model, comparing overlap and staggered fermions. 
They have so far looked at only one lattice spacing, and so far only looked at 
condensate- and topological charge-related quantities. Results from overlap and 
staggered fermions are quantitatively different. However, agreement between stag- 
gered and overlap simulations is considerably improved by smearing the action. 
This is true both for Nf = 1 and 2. There does not seem to be anything funny 
happening with the Nf = 1 fractional root. It may just be that the differences 
between staggered and overlap fermions are cutoff effects, nothing more. The au- 
thors caution that their gauge configurations are much smoother than (my reading) 
typical a — 0.13 fm four-dimensional Nf = 2 + 1 Asqtad ST/ (3) configurations, so 
one should be cautious with interpretations. More work needs to be done (and they 
will do it), varying the lattice spacing and checking scaling. 

All dynamical fermion simulations are costly - at least a factor of 100 more 
than a quenched calculation done at similar values of lattice spacing and quark 
mass. And yet, looking back at Fig. |SJ the change in the numbers which come out 
of the analyses, going from quenched to 2 + 1 flavor QCD, is only a few percent. 
True, these are quantities which are selected to be insensitive to hadronic decay, 
or indeed any strong non quenching effects, but it still seems absurd - to work so 
hard for so little effect. In any classroom physics lecture, this discussion would be an 
introduction to a description of how to do a perturbative calculation, systematically 
improving one's approximate result. But no such story exists for simulations of 

qceJEI 

6. Hadronic Matrix Elements from the Lattice 

One of the major goals of lattice calculations is to provide hadronic matrix ele- 
ments which either test QCD or can be used to provide low energy hadronic matrix 
elements as inputs to test the standard model. A wide variety of matrix elements 
can be computed on the lattice (compare Fig. 112(1 and used to extract CKM mix- 
ing angles from experimental data. A new organization, the Lattice Data Group, 
http : //www . cpt . univ-mrs . f r/ldg/, is attempting to build a "lattice wallet card" 
of relevant results. 

The lattice enters at a rather late stage of the calculation. Consider some 
"generic" electroweak process. Viewed at the shortest distance scale, a distance 
Ax ~ 1/Mwi quarks emit and absorb W and Z bosons. However, physical hadrons 
containing u, d, s, c, or b quarks are big objects, a fraction of a fermi in size, and they 



1, 2008 22:50 WSPC/INSTRUCTION FILE draftl 



Lattice QCD 37 



cannot "see" the highly virtual W and Z bosons. Their interactions are described 
by a low-energy effective field theory valid at scales p, of a few GeV, which would be 
constructed by integrating out short distance physics using the operator product 
expansion, combined with the renormalization group. Th e effective Hamiltonian 
basically reduces to a sum of four-fermion interaction J^S. For example, single- W 
exchange processes would be described by 



eff 
W 



G 

-j=^2ci(ji)Oi(n) 



(75) 



where the GVs are four fermion operators and the c,'s are (known) Wilson coeffi- 
cients. Similarly, the (second-order weak) process B — B mixing is parameterized 
by the ratio x d = (AM) bd /T bd 



Xd 



T bd-^2VQC D F{ 



-)\V; b V td \H(»){-(B\h P (l - J 5 )db lp (l - j 5 )d\B)} (76) 



This is a typical formula used to relate an experimental number to a Standard Model 
prediction. It has a product of factors from phase space integrals or perturbative 
QCD calculations, a combination of CKM matrix elements (whose determination 
is presumably the primary goal), followed by a four quark hadronic matrix ele- 
ment. In order to extract the CKM matrix element from the measurement of x d 
(and its strange partner x s ), we need to know the value of the object in the curly 
brackets, defined as 3/8M bd and parameterized as m 2 B f Bd B bd where B bd is the so- 
called B-parameter, and /b is the B-meson decay constant, (0|&7o75C?|-B) = Jb^b- 
Lattice calculations are the most model-independent way to compute the decay con- 
stants, B-parameters, and their ratio £ = (/b 3 \j Bb s ) / {Jb\J B-q). (Hats denote the 
renormalization group invariant quantities.) The B-parameters and decay constants 
are scale and prescription-dependent quantities and the lattice-regulated quantities 
have to be converted into some continuum regularization scheme to complete the 
calculation. 

The lattice calculation of one of these quantities involves many ingredients. 
First, one has to extract the lattice-regulated matrix element. Most of the matrix 
elements measured on the lattice are extracted from expectation values of local 
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Fig. 12. Matrix elements which can be computed reasonably reliably with lattice methods, and 
their impact on the CKM matrix. From Ref. 61 
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operators J(x) composed of quark and gluon fields. For example, if one wanted 
(0| J(x)\h) one could look at the two-point function 

C JO {t)=^2{0\J{x,t)O(0,0)\0). (77) 

X 

Inserting a complete set of correctly normalized momentum eigenstates 

1 x ^ \A,p)(A,p\ 

ze a (p) (78) 

and using translational invariance and going to large t gives 



(0\J\A)(A\O\0) 



C JO (t) = e- ^< ' - " ■ (79) 
A second calculation of 

Coo(t) = Y,(0\O(x, t)O(0, 0)|0) -> e^'Hf ( 80 ) 

X 

is needed to extract (0| J\A), fitting the two correlators with three parameters. 
Similarly, a matrix element (h\ J\h!) can be gotten from 

C AB (t,t') = Y,(0\O A (t)J(x,1?)OB(0)\0). (81) 

X 

by stretching the source and sink operators O a and Ob far apart on the lattice, 
letting the lattice project out the lightest states, and then measuring and dividing 
out (0\O A \h) and (0\O B \h). 

These lattice matrix elements are not yet the continuum matrix elements. Typi- 
cally, one is interested in some matrix element defined with a particular regulariza- 
tion scheme. It is a generic feature of quantum field theory that an operator defined 
in one scheme (MS) will be a superposition of operators in another scheme (lat- 
tice). In principle, the superposition could be all possible operators. So an operator 
of dimension D will mix like 

</|Or*0«)N>MS =a D £ Z nm (f\O latt (a) m \i). (82) 

The only restrictions on mixing are the ones imposed by symmetries, and generally, 
lattice actions have fewer symmetries than in the continuum, so operator mixing is 
more severe. The most serious source of mixing for light quark operators is the way 
lattice fermions treat chiral symmetry. Wilson-type fermions break chiral symmetry 
(even massless ones do so off-shell) and so nothing prevents mixing into "wrong 
chirality" operators. The mixing structure of overlap fermions (and domain wall 
fermions in the limit of infinite fifth dimension length) is basically identical to the 
continuum, and these formulations really come into their own in cases where chiral 
symmetry is important. 
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In Eq. 1821 the mixing coefficients to lattice operators with the same dimension- 
ality as the continuum ones term will contain the anomalous dimension matrix of 
the continuum operators 

g 2 

Zmn = 1 + -TTT^ilmn log CLfl + A mn ) + . . . (83) 
107T 

This must happen to cancel the mu-dependence of the coefficient function, since 
c(/i)(/|O cont |z) M is independent of the renormalization point. In principle the lead- 
ing log could be summed, but in practice people don't know how much of the 
constant term A should be absorbed into a change of scale of g. This would in- 
volve a two loop calculation. (Brave groupie are beginning such calculations.) So 
they are just left there as constants. There are also terms for mixing with higher 
dimensional operators, which give contributions proportional to positive powers of 
a. (These are usually benign, they look just like scaling violations.) One can also 
have mixing with lower dimensional operators, with contributions involving nega- 
tive powers of a. When they appear, these are dangerous. They must drop out in 
the continuum but it is a delicate business, since they look like they are growing as 
an inverse power of a. 

People compute the ^ nm 's in a number of ways. Most straightforward is to use 
perturbation theory. For lattice actions involving thin links, lattice perturbation 
theory in terms of the bare coupling g(a) is not very convergent: the A nm 's can 
be very large (order 30, sometimes), so a s A/(4ir) ~ 0.2 or more. The source of 
this behavior is the "tadpole graph," usually as part of the fermion self-energy. 
(See Fig. El) The lattice fermion-gauge field interaction looks like l tp(x)U fJl (x)ilj(x + 
ft) and U ~ 1 + igaA^ — g 2 a 2 /2A 2 l + . . .. The ipA'iip vertex, not present in any 
sensible continuum regularization, causes problems when the gluon forms a loop: 
the quadratic divergence from the loop integral combines with the a 2 to give a 
finite contribution-in fact, it is often the dominant contribution. In perturbation 
theory one must also choose the momentum scal e in th e running coupling constant. 
There are reasonable choices for how to do tha ^8 l 73 l Perturbation theory for fat 
link actions is generally better behaved. The fat link vertex acts as a form factor to 
suppress the tadpole. Situations where \Z\ — 1 equals 0.05 or less are not unusual^. 

Often one can find Z nm 's by computing lattice Ward identities"^. For exam- 
ple, with overlap fermions, chiral symmetry "protects" the combination rrtqipip 
and constrains the (local) pseudoscalar and scalar currents, and vector and ax- 
ial vector currents, to have equal renormalization constants. The matrix element 
of the pseudoscalar current gives (O|'075^|PS') = fp S mp s /(2m q ), with no lattice- 
to-continuum renormalization factor for fp S ; fp S = fps- The zeroth component 
of the axial current has as its matrix element (0\'tpjQjc i tp\P S) = fpgirtps, with 
ZaJps — fps- Thus Za — fps/fps which can be computed by fitting a correlator 
with a pseudoscalar sink and a correlator with an axial vector sink, and taking 
the ratio. (To experts, this is slightly wrong: m q r ip(l + aD)ip and the point split 
currents f (1 + aD)ip are the quantities related by overlap current algebra. Experts 
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Fig. 13. The "tadpole graph" for the fermion self energy.. 

also know that the u aD" terms are removed by a slight redefinition of the lattice 
propagator...) 

The most commonly used nonperturbative method for comp utin g matching 
factors is called the "Regularization Independent" or RI schemed One com- 
putes quark and gluon Green's functions in a smooth gauge, regulated by giv- 
ing all external lines a large Euclidean squared momentum p 2 . For example, for 
a bilinear, one might compute an un-amputated correlation function Go{p-a). 
One would also compute the quark propagator in momentum space, 5(p, a). 
Then the vertex would be found by snipping the propagators off the correlator, 
Ao(p, a) = 5 _1 (p,a)C7o(p, a)5 _1 (p,a) and the vertex would be computed by pro- 
jection: To (p, a) (x TtAoPo ■ Roughly speaking, the Z factor is the ratio of To {p, a) 
to its free-field value, or one would match to a perturbative calculation of the ver- 
tex, regulated identically. To measure a Z-factor, one should see a plateau in the 
observable versus p 2 . One must stay away from large pa (where discretization effects 
enter) and small p (where nonperturbative physics turns on), The author has not 
tried doing this himself, but there are many claims in the literature from people 
who have, that it is straightforward to do. 

People become very dogmatic about how to compute matching factors. It is a 
very interesting exe rcis e to compare the points of view of a firm advocate of lattice 
perturbation theorjE^I with that of an equally firm advocate of nonperturbative 
renormalizatiorP^. Read these papers back to back! 

6.1. Charm and bottom matrix elements 

We spent a long section on relativistic fermions, so the reader will not be surprised 
to lear n th at three are also many choices for discrctizing heavy fermions. It is very 
difficult^ to do direct calculations with relativistic b quarks because the lattice 
spacing is much greater than the b quark's Compton wavelength (or the UV cutoff 
is below rrib). However, it might be better to think of the lattice theory as an effective 
field theory for low-momentum physics in the presence of two high energy scales-the 
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cutoff (inverse lattice spacing) and the heavy quark mass. The effects of the short 
distance are lumped into coefficients of the effective theory ^ . A simple example 
of such a theory would be the standard clover action, since it has lattice versions of 
every operator with dimension < 5. So the simulation involves a heavy clover quark 
and (perhaps) some other type of light quark. The heavy quark dispersion relation 
is E(p) ~ mi + p 2 /(2m 2 ) + . . . (if p << m) and it might happen that mi 7^ m 2 
because of discretization artifacts. The heavy quark magnetic moment is fi = l/m^ 
and again, m 3 ^ m 2 ^ m\. But mi is just an overall constant, so one could shift the 
measured hadron mass by m 2 — mi and tune the shifted mass to the desired value 
(to do the simulation at the physical B mass, for example). The lattice magnetic 
moment could also be tuned (by varying the coefficient of the clover term) using 
(say) the B* — B hyperfme splitting. Hadronic parameters not used to determine 
the lattice parameters are then objects of calculation. 

Nonrelativistic QCD has also been discretized and used to make very precise 
calculations of the properties of quarkonia^l This formalism can also be used for 
the heavy quark in a B or D meson (again as long as its momentum is small.) The 
"static" limit (infinite 6-quark mass) is often used as an additional point on the 
curve. Then one can try to interpolate all the way from light quarks to heavies and 
get all the decay constants at once. 

Another new de velop ment is the use of staggered quarks for the light quark in 
a heavy-light systeirP-4 This will allow one to get to smaller quark masses than 
Wilson- type valence quarks permit, which will be important for chiral extrapola- 
tions. Overlap light quarks are even prettieJSSI but much more costly. Calculations 
using overlap fermions will probably appear next yea r. 

Early form of fat links failed rather spectacularljE^^I when used with heavy 
flavor simulations because they smeared away the short distance part of gluonic 
interactions, but HYP links have been used to enhance the signal of static-light 
correlator:^!. 

An important systematic effect is the use of unphysical (too heavy) values of 
the light quark masses both in quenched and unquenched simulations. People had 
just been doing linear extrapolations to get down to the physical light quark mass, 
but this is not correctr ^*^"^! One has to use chiral perturbation theory, with its 
non-analytic (logarithmic) behavior at small m q . An example of such behavior is 

^-l = (m%-ml)f 2 (A) 

JB 



l + 3g 2 



1 1 3 

-I(m K ) + -I(m v ) - -I(m n ) 



(84) 



wherelp (mps) — rn? PS ln(mp S /A 2 ) and / 2 is a low-energy constant. Kronfeld and 
RyaifSl pointed out that the inclusion of chiral logarithms in the chiral extrapola- 
tion of lattice data for heavy-light dec ay c onstants can drastically change /b b //b 
and £. By assuming g 2 = g 2 D , DlT — 0.35^1 and / 2 = 0.5(3) GeV -2 , they concluded 
that £ = 1.32±0.10, which is more than 10 per cent larger than Ryan's 2001 global 
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lattice estimate^! 

Most lattice results are done in quenched approximation. In many cases, all the 
other systcmatics (extrapolation to the continuum, matching factors, etc.) can be 
beaten down to below a few per cent. But the quenching systematic won't go away 
until dynamical fermion simulations are done. 

There are many lattice calculations of fg, fjj, Bb , and form factors for semilep- 
tonic decay. My summary leans rather heavily on reviews by Wittig*^, BecirevkP^I 
and KronfelcP^I. 

The decay constant is computed by combining a heavy quark and a light anti- 
quark propagator into Eq. 1771 Decay constants probe very simple properties of the 
wave function: in the nonrelativistic quark model /m = ' t P(^)/^/ m M, where ip{0) is 
the qq wave function at the origin. For a heavy quark (Q) light quark (q) system 
ip(0) should become independent of the heavy quark's mass as the Q mass goes 
to infinity, and in that limit one can show in QCD that ^ttlmSm approaches a 
constant. 

The lattice calculation of decay constants which I know best is by the MILC 
collaboration, again, Refs. (Nf = 2) and^O (Nf = 2 + 1, preliminary results). 
These authors have done careful simulations, in quenched, 2-flavor, and 2+1 flavor 
QCD, at many values of the lattice spacing, which allows one to extrapolate to the 
continuum limit by brute force. Examples of results are shown in Figs. Upl and 1151 
The first figure shows fVM vs 1/M, with data extending from the static limit to 
below the D. One would fit this data to a power series in 1/M to get to the B or 
D. Notice that the simple heavy quark assumption ju oc l/VM fails by 50 per 
cent at the B. 

A plot of the data which goes into Fig. El is shown in Fig. El There does 
not seem to be a lot of lattice spacing dependence. Quenched data extends to 
smaller lattice spacing, and so one might want to use its slope to extrapolate to 
the continuum. The two a — points show different extrapolation choices, clearly 
a source of uncertainty. 

Figure HHI from Ref. llOOl shows how fs changes as the sea quarks switch on - a 
jump from 170 to 210 MeV. The ratio /b s //b, shown in Fig. 1171 seems to be much 
more stable. Sometimes, apparently, quenching works well, sometimes not. 

Semileptonic decays involve evaluating lattice correlators of three-point func- 
tions (two mesons and the current). One just fits the data (measured hopefully at 
many masses and as many momenta as possible-remember, momentum is quantized 
to be a multiple of 2tt/L in a box of size L) to the expected set of form factors. For 
example, in B — > niv, 

(n(p)\V,\B(p>)) = f + (q 2 )[p'+p- m% ~ ml q], + f{q 2 ) ml q2 ml q ,. (85) 

The best signals come when the momenta of the initial and final hadron are small. 
Then the large B mass forces q 2 (q = lepton 4-momentum = ps — Pn) to be large. If 
the form factor is needed at q 2 ~ 0, a large extrapolation is needed, and there will 
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Fig. 14. Pseudoscalar meson decay constant vs 1/M, from Ref. 
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Fig. 15. Jb vs. a from Ref. 90 showing extrapolations to the continuum limit of full QCD data. 



be additional errors and model dependence in the answer. (Lattice people have no 
advantage over anyone else at guessing functional forms.) However, finding V u b from 
experimental data only requires knowing the form factor at one value of q 2 . This 
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Fig. 16. /b as a function of dynamical quark mass (plotted as (m^ /m p ) 2 ) from Rcf. 100 There 
appears to be a considerable shift as the sea quarks drop in mass. 
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Fig. 17. /s.//s as a function of dynamical quark mass (plotted as (m^/m p ) 2 ) from Rcf. 100 
This ratio is almost independent of the sea quark mass. 

should work so long as the experiment has enough data to measure the differential 
rate around that region of q 2 . 
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As an example of a recent approach, the FNA L group has measured B — ► 
Dlv and B — > D*lv form factors at zero recorP^. They have a clever technique 
for removing much of the lattice-to-continuum Z-factors by computing ratios of 
matrix elements. For example, in B — > D£v, the differential cross section at zero 
recoil invo lves a form factor h +(v • v' = 1). This quantity is one in the heavy 
quark limilPSl with correction :]-*-^ ! which start at X/m\. The lattice calculation 
uses lattice operators and needs to be matched to a continuum result. The FNAL 
group extracted it from the dimensionless ratio 

(D\c lQ b\B)(B\b l0 c\D) , 2 

— — - - - = \h+{v -v = 1) (86) 

{D\c-/ c\D){B\tr/ b\B} 

The denominators are just diagonal matrix elements of the charge density, (i.e. they 
are basically just charges) and the ratio of renormalizations between numerator 
and denominator is under better control than the individual factors. The B — > 
D*£v calculation is similar. Their calculation of its form factor is h(l) — 0.913 
with a combined non-quenching error of about two per cent. Using this number 
in conjunction with experimental measurements of the branching ratio would give 
about a five per cent uncertainty in \V c b\, not including the quenching uncertainty, 
of course. 

Finally, £. I know of two published calculation of £ in Nf = 2 QCD, 
which inclu des c hiral logarithms in their analysis. The first is by the JLQCD 
collaboratioiffSI. Their quark masses (0.7m s < m q < 2.9m s ) are too large to 
actually see any chiral logarithm, but they can boun d the size of the effect through 
their inability to see it. Recently, Wingate et aP^ presented results (from simula- 
tions using Asqtad valence quarks on the MILC Asqtad background configurations) 
which are consistent with the behavior expected from chiral logarithms. A value 
of £ well away from unity is suggested, about 1.25. Kronfeld's summary is shown 
in Fig. [TBI The light quark data of Ref. EH tracks the fit of Ref. EH However, 
it happens that the decay constants of the two groups are quite different: about 
260 MeV for Ref. ESI as opposed to 215 MeV for Ref. EH Errors in both cases 
are claimed to be about ten per cent from all sources. Clearly something is not 
working... A glance at Fig. 1161 shows a third prediction for Jb of about 210 MeV 
(from Jb, see Fig. 11611 . xl.2 for £ from Fig. El — 250 MeV, with maybe a ten per 
cent error. 



6.2. Kaon Matrix Elements 
6.2.1. B K 

The kaon B-parameter Bk 1 defined as the matrix element of a particular four 
fcrmion operator + , 

l(m K f K ) 2 B K = (K\s^(l - 75)^(1 - 7 5 )d\K), (87) 
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Fig. 18. £ vs quark mass. The black points are from Ref. 104 while the grey ones are from Ref. 
[Tulfl The straight line is a linear extrapolation and the dotted line is a chiral logarithm fit, both 
from Ref. I1U4I This picture is from the review of KronfelcHH 



is an important i ngred ient in the testing of the unitarity of the Cabibbo-Kobayashi- 
Maskawa matrix ' ^" 1 It has been a target of lattice calculations since the earliest 
days of numerical simulations of QCD. Bk puts a severe test on the chiral behavior 
of lattice simulations and so it is a useful test of methodology (in addition to being 
a dimensionless number interesting to experiment). There has been a continuous 
cycle of lattice calculations using fermions with ever better chiral properties. Un- 
fortunately, to date all calculations (except one preliminary one described below) 
have been done in quenched approximation, and so there is no lattice prediction 
relevant to experiment. 

The first lattice calculations of Bk were done with Wilson-type actions. These 
actions break chiral symmetry. This is a problem, because while the matrix element 
of + scales as Wpg, the opposite chirality operators which appear under mixing 
have matrix elements which are independent of mps- The signal disappears under 
the background in the chiral limit. Techniques for handling operator mixing have 
improved over the years, (for recent results, see Refs. 110711 but I think this 
approach is still very difficult. 

Staggered fermions (Refs. Ill 1H112|) have enough chiral symmetry at nonzero 
lattice spacing, that operator mixing is not a problem. One can obtain extremely 
precise values for lattice-regulated Bk at any fixed lattice spacing. However, to 
date, all calculations of Bk done with staggered fermions use "unimproved" (thin 
link, nearest-neighbor-only interactions), and scaling violations are seen to be large. 
For example, the JLQCD collaboratio n I saw a thirty per cent varia tion in Bk 
over their range of lattice spacings. Fat link calculations are in progresJU^l 

For domain wall fermions, chiral symmetry remains approximate, though much 
improved in practice compared to Wilson- type fermions. Two groupJU^ESI have 
presented results for Bk with domain wall fermions. Ref. II 141 has data at two lattice 
spacings and sees only small scaling violations. There is a few standard deviation 
disagreement between the published results of the two groups. 
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Fig. 19. Bk comparisons vs lattice spacing, from a selection of simulations with reasonably 
small error bars. Results are labeled diamonds and fancy diamond, Ref. 11121 the fancy cross, Ref. 
11171 octagons, Ref llX4| the cross, Ref. 111^1 and squares, Ref. Illjjl 



Figure EH shows results for Bk at various lattice spacings, for a selection of 
simulations which have reasonable statistics and small error bars. Figure l^Ul prcscnts 
results which are either extrapolated to the continuum limit, or presented by their 
authors as having small lattice spacing artifacts. 

Bk shows quite a bit of nonlinear quark mass dependence (compare Fig. l21l for 
a typical recent result). This suggests strongly that chiral logarithms are present 
in the data, and it is reasonably easy to fit the data including their effects, out to 
quark masses which are about equal to the strange quark's mass. It happens that 
the coefficient of the chiral logarithm is the same in full and quenched QCD. 

The only group reporting results for a calculation of Bk in unquenchedQCD 
that I am aware of is the RBC collaboration, using domain wall fermion J I At 
this year's lattice conference, they presented a preliminary number of B^ DR (/i = 
2 GeV) = 0.50(2) for Nf = 2. It is low compared to all the quenched Bk results 
quoted above. This is interesting: a too-low value of Bk moves the region in the 
(p, iff) plane of the CKM matrix allowed by Bk (actually, ck ) away from the heavy 
flavor determinations. The authors do not claim anything so grand; their calculation 
is still in its early stages. Phenomenologists should keep an eye on their result, 
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Fig. 20. Bk comparisons presented "as if" they were taken to the continuum limit. The label 
in parentheses characterizes the kind of lattice fermions used: W for Wilson, CI for Clover, KS 
for staggered, DW for domain wall, and OV for overlap fermions. References are CBS, Ref. Iltl/I 
JLQCD(W), Ref. [TT0] Conti, Ref . [MD KGS, Ref . [TTT] JLQCD(KS), Ref . [TT2] CP-PACS, Ref. 

rrm rbc, Ref. rmn gghlr, Ref. rrm and milc ittsi The points of Rets, rrrni UTTi rrrei [m 

andlTTs! are the results of a a continuum extrapolation; all the rest are simulations at one lattice 
spacing. 



however. 



6.2.2. K -> 7T7T and e'/e 

There are two types of CP violating decays of kaons. The first (and largest) origin 
is the admixture of a small amount of a CP eigenstate into the wave function of 
a particle which is dominantly of the other eigenstate. This is called "indirect CP 
violation," as opposed to 'direct CP violation," in which the decay of the large 
component of the wave function violates CP. The ratio of these two processes is 
parameterized by the ratio 

- = (16.6 ± 1.6) x 10~ 4 (88) 
where I am quoting its 2003 experimental value, from the review of BuraJEO] 

based 

on experimental data of Refs. IT^ rEffl IT231 

About two years ago, the RB<ji±£l and CP-PACS 126 collaborations computed 
e'/e in quenched approximation, using domain wall fermions. Their numbers are 
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Fig. 21. A typical (not very good) lattice calculatior P^ lof B^ DR (fi = 2 GeV) as a function of 
quark mass from the two data sets-diamonds for a = 0.13 fm, squares for a = 0.09 fm. 



in spectacular disagreement with Ea. 1881 Ref. 11151 quotes —7.7(2) x 10~ 4 and Ref. 
126 says that their value is "negative, and has a magnitude of order 10~ 4 . What is 
going on? 

Lattice calculations of e'/e are complicated. They begin with the low energy 
effective Lagrangian of Ea. 1751 to compute (first) K — > irir amplitudes 

K° -> (tttt)/ = Aiexp(iSi) (89) 

and the parameter related to the AI =1/2 rule, 

1 - ^444 - 22. (90) 

These feed into the expression for e'/e, which happens to be dominated by the 
difference of two of the operators, the QCD penguin 

Qe = (sd) L ((uu) R + (dd) R + (ss) R ) , (91) 

and the electroweak penguin 

Qs = \{sd) L ({2uu) R - (dd) R - {ss) R ) , (92) 
-MVuKt) x (93) 



e 2|e|Re^ 

2/6^)(Q 6 )tS(M) - ^8(M)(Q 8 )]g(M) ) (i)-ll 
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and {Qi)iexp(i5i) = {iriri\Qi\K). 

People don't calculate K — > nn directly on the lattice; it is difficult " to extract 
the phase shifts from the mr final state interactions from lattice data (never mind 
trying to separate the two pio ns to asymptotically great distances). Instead, they 
use chiral perturbation theory ESI to relate K — > tttt amplitudes to K — > ir and 
K — ► vacuum. In the case of the AI = 3/2 amplitude there is a factor of two change 
in the lattice result depending on whether tree level or one loop chiral perturbation 
theory is used. This is shown in an old figure from Ref. |128 in Fig. 1221 
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Fig. 22. AI = 3/2 K — > 7™ amplitude with (fancy symbols) and without (plain symbols) one- 
loop corrections of quenched chiral perturbation theory. Data are crosses and fancy crosses, Ref. 
11301 diamonds and fancy diamonds, Ref. ll.3l1 



At the end of the day, there is also a conversion from the lattice regulated matrix 
element to the continuum-regulated ones. This can be done either perturbatively 
(worrying about how convergent the calculation is) or nonperturbatively. 

The quenched approximation is a real problem. Recall that the flavor symmetry 
of quenched QCD is not SU(3) L <g) SU(3) R , it is the graded algebra SU(N\ N) L ® 
SU(N\N)r. Consequently, the flavor content of penguin operator Oq is altered: it 
came from running short distance physics down to long distance using unqucnchcd 
QCD evolution, but then it is going to be evaluated in quenched approxi matio n. 
Extra low energy constants appear, which would not be present in full QCEff=2l 
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As a practical consequence, matrix elements of these operators have a chiral 
extrapolation similar to that of Bk , 



£m 2 



< TrIO, | A' >= C(l + H™%s)) + bm 2 PS . (95) 



Golterman and Pallentff^ computed the coefficients for chiral logarithms in 
quenched and partially-quenched QCD, and found that £ = for Os (in quenched 
approximation and in the degenerate-mass limit). This is different from in full QCD. 
After seeing the rather large curvature of Bk in Fig. 1211 plausibly effect of a chiral 
logarithm, this discovery is cause for unease. 

Another peculiarity of the quenched approximation is that the B-parameters for 
many of these operators go to zero in the chiral limit, because mp s /m q diverges. 
Recall that the B-parameters for these operators are defined as the ratio of the 
operator to its vacuum-saturated approximation. For left-right operators, 

B 3/2 = (A|Op|7r)_ (96) 

c P (it I V>75 i> 1 0) (0 1 ip"f5tp | K) + ca (tt I ^7^75^ 1 0) (0 1 V'Tm 75 tp I A") 

where the Cj's are numerical coefficients. The PC AC relation says that 
(0\ipj5ip\P S) = \rn 2 psl m qfps- The divergence of (m ps /m q ) 2 in the denomina- 
tor is not compensated by any singular behavior in the numerator. 

My conclusion is that we should be happy that people are building all the tools 
to do these calculations, and that one shouldn't take any of the quenched numbers 
seriously. RBC and CP-PACS are doing simulations in full QCD, and it will be 
interesting to see what they find when they revisit e'/e in the next few years. 
I suspect that e'/e will teach us more about how to do calculations than about 
whether the Standard Model is right or not. 



7. Instead of conclusions 

Lattice QCD is in a period of rapid evolution. The quenched approximation is 
fading out as a vehicle for precise calculations of quantities to be directly compared 
with experiment. It was helped on its way by algorithms for simulating QCD while 
preserving flavor and chiral symmetry, which cause one to attempt to do problems 
whose analysis forces one to confront the fundamental differences between quenched 
and full QCD. Simulations are beginning to reveal the presence of quenched chiral 
logarithms at small quark mass. 

Full QCD simulations are being used in phenomenology as never before. The 
authors of Ref. 1611 are carrying out an extensive program of light and heavy fla- 
vor physics using Asqtad dynamical fermions. They argue that essentially all the 
entries of the CKM matrix (except Vtb) can be impacted by lattice calculations 
whose systematics are well understood. Asqtad fermions will eliminate the unknown 
quenching systematic. Results from these simulations are already appearing, and 
more will appear within a year. 
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A big part of the lattice QCD computing resources in this country will go into 
the construction of large sets of Asqtad lattices, which the idea of making the 
data sets publicly available. Other people will use their configurations as "back- 
ground fields" for other matrix element calculations. In many cases the cost of the 
simulations which are needed directly to compute a particular matrix element is 
much smaller than the cost of generating the configuration with its dynamical light 
quarks. Quenching systematics would presumably be a thing of the past. These 
calculations promise fantastic improvement in the quality of lattice determinations 
of hadronic matrix elements. 

The MILC collaboration is also planning to push to smaller lattice spacing and 
to smaller nonstrange quark mass. Recall that their present simulations are done at 
Tn Ut d/Tn s — 1/5, while in the real world the ratio about 1/20. The smallest lattice 
spacing in their simulations is about 0.09 fm. They estimate that the cost of this 
data set is about 0.8 Petaflop-hours. The cost of dropping the quark mass by a 
factor of two at fixed physical volume, or of dropping the lattice spacing at fixed 
quark mass, as given by Eq. 1301 are each less than a factor of 16 increase over 
present simulations. They believe that a combination of these two simulations will 
allow extrapolations to the continuum and chiral limit at the one per cent level. 

This won't happen in a year, but it is not completely unreasonable to think 
about doing it. The computer resources will be there in the next decade. 

However, I think that without some better understanding of the "det 1 ^ 4 " re- 
placement, all these results will be clouded by uncertainty: are these simulations 
fundamentally correct, or not? And I don't think that the answer to this question 
will come from numerics, at least not "direct" numerics in four dimensions, like 
spectroscopy or matrix element calculations. I believe it is an issue of principle. 

Dynamical simulations with domain wall or overlap fermions need a lot of de- 
velopment before they can be attempted with finite computer resources. 

Lattice QCD at the end of 2004 promises to be quite different from lattice QCD 
at the end of 2003. 
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